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ABSTRACT 

We present the broadest and most precise empirical investigation of red giant 
branch bump (RGBB) brightness and number counts ever conducted. We imple- 
ment a new method and use data from two Hubble Space Telescope (HST) globu- 
lar cluster (GC) surveys to measure the brightness and star counts of the RGBB 
in 72 GCs. The measurement precision reaches 0.005 mag in brightness and 10% 
in the number counts. The position of the main-sequence turnoff (MSTO) and 
the number of horizontal branch (HB) stars are used as comparisons where ap- 
propriate. Several independent scientific conclusions are newly possible with our 
parametrization of the RGBB. Both brightness and number counts are shown to 
have second parameters in addition to their strong dependence on metallicity. 
The RGBBs are found to be anomalous in the GCs NGC 2808, 5286, 6388 and 
6441, likely due to the presence of multiple populations. Finally, we use our 
empirical calibration to predict the properties of the Galactic bulge RGBB if the 
assumption of similar stellar physics for the bulge and Galactic GC system holds. 
The updated RGBB properties for the bulge are shown to differ from those of 
the Galactic GC system, with the former having lower number counts, a lower 
brightness dispersion and a brighter peak luminosity than would be expected 
from the latter. This discrepancy is well explained by the Galactic bulge having 
a higher helium abundance than expected from GCs, AY~ +0.06 at the median 
metallicity. 

Subject headings: Hertzspring- Russell and C-M diagrams - Galaxy: globular clus- 
ters: general - Galaxy: globular clusters: individual: (NGC 2808, NGC 5286, 
NGC 6388, NGC 6441) - Galaxy: Bulge 
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Introduction 



The red giant branch bump (RGBB) is a prominent feature of color-magnitude diagrams 



(CMD) along the red giant (RG) branch that was first theoretically described by Thomas 



(1967) and Iben (1968). It was first empirically confirmed nearly two decades later, by King 
et al. (1985), in their observations of the Galactic globular cluster (GC) 47 Tuc. During the 
first ascent of the RG branch when the hydrogen-burning shell moves outward, stars become 
temporarily fainter before becoming brighter again due to a discontinuity in the chemical 



abundance profile near the convective envelope (Sweigart et al. 1990). As the star thus 
has the same luminosity on three separate occasions, an excess in the luminosity function 
at a characteristic magnitude becomes visible in the luminosity function of RG stars. The 
properties of this excess, such as its characteristic brightness and expected number counts, 
are a steeply sensitive function of a stellar population's age, initial helium abundance, and 



metallicity (Cassisi & Salaris 1997; Bono et al. 2001; Bjork & Chaboyer 2006; Di Cecco et 



al. 


2010 


Cassisi et al. 


2011; 


Nataf et al. 


2011a|b 



Even though the RGBB occurs before more complex phases of stellar evolution such 
as the helium flash, and should as such be well-constrained theoretically, there has been an 
ongoing debate in the literature as to a discrepancy between the predicted and observed 



brightnesses of the RGBB. In their pioneering study of the RGBB in 11 GCs, Fusi Pecci et 



al. ( 1990 ) found that a zero-point shift of several tenths of a magnitude was required to bring 



theory into agreement with observations - the predicted luminosity of the RGBB was greater 



than that observed. Some disagreement followed. Cassisi &; Salaris (1997) and Zoccali et al. 



(1999) argued that there was no significant disagreement between theory and observations if 
one accounted for observational uncertainties in [a/Fe] and [Fe/H] abundances. Moreover, 
as the the brightness of RGBB was measured relative to the HB, it was not clear how any 
discrepancy between models and data should be interpreted. 
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Fig. 1. — TOP: The values of My^GBB and their statistical errors for all 72 Galactic GCs 
studied in this work are shown as filled red circles. The error bars are the statistical errors 
due to the fit, and do not include the systematic error from the assumed extinction and 
distance values. RGBB properties in the anomalous GCs NGC 5286, 2808, 6388 and 6441 
delineated by black x-shaped symbols. BOTTOM: Number counts on the RGBB normalized 
with respect to the number of HB stars as a function of [M/H]. The values of Jjigbb f° r the 



Galactic spheroid systems M32 and the MW bulge, respectively measured by Monachesi et 



al. (2011 ) and Nataf et al. (2011a), are shown as the empty black squares in the lower panel. 



However, since those works the breadth of photometric catalogs, the range of available 
diagnostic, and the accuracy of metallicity scales have evolved. More recently, three works 
using three distinct methods found that the brightness of the RGBB is 0.2-0.4 mag fainter 
than predicted by models. Di Cecco et al. (2010) used a large sample of 15 GCs spanning 
2 dex in metallicity. They found that the difference in l/-band between the RGBB and the 
zero-age HB, ISV^ahb = Vrgbb — Vzahb, was larger than predicted by models, i.e. that the 
bump was fainter in observations than in models. The discrepancy was found to be much 
more severe in metal-poor GCs. Cassisi et al. (2011 ) compared the RGBB brightness relative 
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to the main sequence turnoff (MSTO) in the Hubble Space Telescope (HST) F606W band, 
AF606W£ggj? = F606W MS to - F606W RGBB . They found that agreement with theory 
could only be obtained if GCs are in fact ~4 Gyrs younger than the ages resulting from 



distances inferred from main-sequence fitting. Troisi et al. (2011) investigated the empirical 



brightness difference between the RGBB and the point on the main sequence that is at the 
same color as the RGBB, which they label A£. Models predict this to be a very sensitive 
helium diagnostic that only weakly depends on age. Unfortunately, these same models only 
match the data if the helium abundance in the clusters is set to an unphysical value of 
Y=0.20 - lower than the primordial value Y=0.249 derived from observations of the cosmic 



microwave background and the standard model of particle physics (Steigman 2010). Three 



different methods lead to one consistent result, that stellar models overestimate the predicted 
luminosity of the RGBB. 

A second observational challenge in RGBB astrophysics has also recently come on the 
scene. Nataf et al. (2011a) found that the RGBB in the Galactic bulge had deficient number 



counts relative to both stellar models and observations of GCs. They measured an RGBB 
to red clump (RC) ratio f^GBB = 0.13 ±0.02 toward Baade's window, compared to ~25% in 
the most metal-rich GCs and stellar model predictions that the lifetime of the RGBB at the 
age and metallicity of the bulge should be ~25 Myr (Nataf et al. 2011a), which predicts a 
~25% number fraction if one assumes a lifetime of ~100 Myr for the HB (Han et al. 2009). 



Meanwhile, in an analysis of HST CMDs toward the dwarf elliptical M32, Monachesi et al. 
(2011) measured Jjigbb = 0.154±0.036. This measurement was done toward a similar stellar 
system, but with completely distinct photometry, systematics, and a different fitting method 



than that used in Nataf et al. (2011a). Bono et al. (2001) and Di Cecco et al. (2010) both 



mentioned that enhanced helium would reduce the lifetime of the RGBB, but the discussion 
was not broad nor deep enough to quantify how much helium would be necessary, and if 
a consistent solution was available to bring in line each of brightness, brightness dispersion 
and star counts. Further investigation was required to precisely ascertain the nature of this 
unexpected finding. 

Hence the epistemological basis for this paper. We wish to empirically calibrate the 
RGBB as a probe of the nature of the stellar populations from which they emerge - com- 
parison to stellar evolution models will follow in a subsequent paper (Nataf et al. 2012). 
We argue that the RGBB is one of the most undertapped features of the CMD with which 
to interpret stellar populations, and that its diagnostic power will only grow in the next 
decade. We present both an optimal strategy to measure the parameters of the RGBB, and 
an investigation of how the RGBB behaves in the 72 GCs probed by this work. To that 
end, we state our key empirical results on the brightness and number counts of the RGBB 
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in Galactic GC directly in the introduction: 

M V!RGBB = (0.600 ± 0.013) + (0.737 ± 0.024)([M/H] + 1.110), 5 = 0.077 (1) 

a Vrgbb = Vmsto ~ V RGBB = (3.565 ± 0.012) + (-0.549 ± 0.023)([M/H] + 1.152), S = 0.072 

(2) 

Alf$jj B = Irgbb ~ Irhb = (0.123 ± 0.018) + (0.852 ± 0.045) ([M/H] + 0.640), 5 = 0.051 

(3) 

EW RGB B = (0.248 ± 0.010) + (0.121 ± 0.018)([M/H] + 1.134), S = 0.000 

Sbgbb = N ^ BB = (0.111 ± 0.005) + (0.109 ± 0.011)([M/H] + 1.273), 5 = 0.018, (4) 

Nhb 

where the five parameters denote the absolute magnitude of the RGBB peak in V, the bright- 
ness difference in V between the RGBB peak and the position of the MSTO, the brightness 
difference in / between the RGBB peak and the mean brightness of red horizontal branch 
(RHB) stars, the equivalent width EWrgbb of the RGBB feature in the RG luminosity 
function relative to the underlying RG continuum, and the relative number of RGBB stars 
to HB stars. The values of 5 refer to the intrinsic scatter in these relations that is in addi- 
tion to the scatter due to statistical errors in the measurements and the errors in the input 
spectroscopic metallicities. The significant values of S demonstrate that all but one of these 
variables has at least one second parameter. The problematic GCs NGC 2808, 5286, 6388 
and 6441 are left out of the fit. The two key empirical trends probed by this work, quantita- 
tively summarized in the above equations and graphically summarized in Figure [TJ are the 
decreasing luminosity and increasing number counts of the RGBB as metallicity is increased. 

In this paper, we summarize the input photometric and spectroscopic data as well as 
necessary calibrations used in Section [2] A fitting procedure for the different observational 
parameters of the RGBB including star counts is outlined in Section [3j The methodology 
with which we measure the relevant parameters for the MSTO and the HB are discussed 
in Sections [4] and [5] Results for the brightness, number counts and other parameters are 
respectively discussed in Sections [6j [7], [8j A more detailed analysis of the anomalous RGBBs 
in the clusters NGC 2808, 5286, 6388 and 6441 is to be found in Section [9j We use the 
metallicity distribution of bulge stars to derive what the properties of the Galactic bulge 
RGBB toward two distinct sightlines would be if bulge stars have the same input physics 



as Galactic GCs stars in Section 11 Our conclusions are presented in Section 12 We use 
Monte Carlo methods to demonstrate the reliability of our approach in the Appendix. Tables 
summarizing the measurements are to be found following the Appendix. 
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Data 



In this study we make use of two different large-scale GC surveys conducted with HST. 
We also use ground-based data for comparisons to the Galactic bulge. 



We use photometry obtained with HST's Advanced Camera for Surveys (ACS) (Saraje 



dini et al. 2007 Dotter et al. 2011), hereafter "the ACS GCs". The data were taken as part 



of an HST treasury program to obtain high signal-to-noise ratio photometry down to the 
lower main sequence for a large number of Galactic GCs. Artificial star tests demonstrate 
that the photometry is expected to be very precise and complete at the brightness of the 



RGBB (Anderson et al. 2008). 



We also use the HST F439W and F555W photometry obtained as part of an HST 



GC survey program with the WFPC2 camera flPiotto et al.||2002[ ), hereafter "the WFPC2 
GCs". We use the photometry from this database for 16 clusters with well-populated red 
giant branches that were not observed within the ACS survey. Combining these two datasets 



yields a richer sample with better completeness over the metallicity range of GCs. Zoccali 



et al. (1999) and Riello et al. (2003) have already studied the RGBB in the WFPC2 GC 



sample, however we wish to study the two samples together, using a uniform methodology. 



The metallicities for these clusters are taken from the GC metallicity scale of Carretta 



et al. (2009b), except for the three GCs studied in this work not listed there. The metallicity 



for Lynga 7 is taken from Bonatto & Bica (2008), and those of NGC 6426 and Pyxis are 



taken from Dotter et al. (2011). The remaining GC parameters are taken from the Harris 



(1996, 2010 edition) catalog of parameters for Milky Way GCs. 



OGLE-III optical photometry is used to compare the RGBB in GCs to that in the 
Galactic bulge. OGLE-III observations were obtained from the 1.3 meter Warsaw Telescope, 
located at the Las Campanas Observatory in Chile, and are complete to magnitude ~20.5 
in both V and /. Detailed descriptions of the instrumentation, photometric reductions and 



astrometric calibrations are available in Szymahski et al. (2011) 



2.1. Some Calibrations to the Input Data 

We make two adjustments to the input data motivated by the need for uniformity. The 
first is to the definition of V used in the WFPC2 dataset, and the second relates to the 
metallicity scale. 

With respect to the photometric calibration, neither dataset obtained data in V. The 
WFPC2 dataset obtained photometry in F439W and F555W, whereas the ACS dataset 
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photometry is in F606W and F81AW. Photometric values respectively transformed into the 
(V, B — V) and (/, V — I) plane were given. For 13 of the WFPC2 clusters that had also 
had data in the ACS dataset, namely NGC 104, 362, 1851, 2808, 5904, 5927, 5986, 6304, 
6388, 6441, 6624, 6637, 7078, we compared photometric values obtained for the two datasets 
at the level of the RGBB. We found that the V given in the WFPC2 dataset was ~0.0365 
mag fainter than the V given in the ACS dataset. We adjusted the definition of V in the 
WFPC2 results, without adjusting the definition of (B — V). 



The second calibration pertains to the metallicity. The metallicity scale of Carretta et 



al. (2009b) has values of [Fe/H] and ope/H] f° r every cluster studied in this work. However, 
it does not have values of [a/Fe] for all the clusters. We computed a linear fit to all clusters 
in the catalog that have both [Fe/H] and [a/Fe] values, and are not associated to dwarf 
spheroidal galaxies (Arp 2, NGC 4147, NGC 6715, Palomar 12, Terzan 7, and Terzan 8). 
We obtained [a/Fe] = 0.342 — 0.033([Fe/H]+l). This was the value used for all clusters that 
do not have an [a/Fe] value. 

Measurement errors in [Mg/Fe] and [Si/Fe] were derived for 17 GCs spectroscopically 
investigated in Carretta et al. ( 2009a[ ). We take the mean of these errors as the estimate for 



C[a/Fe]- The mean error is <J[ a /Fe] =0.060. Since the scatter in the assumed relation of [a/Fe] 
= 0.342 - 0.033([Fe/H]+l) is 0.085 dex, we assume an error of a [a/Fc] = V0.085 2 - 0.060 2 = 
0.060 dex for all the remaining cluster for which we don't have a reported value of [a/Fe], 
and thus no measurement error. It is likely a coincidence that the mean measurement error 
in [a/Fe] is equal to the intrinsic scatter to the [Fe/H]-[a/Fe]relation. 



We then compute the "total metallicity" of the GCs using the relation of Salaris et al. 



(1993): 



[M/H] = [Fe/H] + log(0.638 * 10 [a/Fc] + 0.362) (5) 
This relation is equivalent to the statement: 

10 [M/H] « 0.3Q2(N Fe /N FeQ ) + O.Q38(N a /N a(D ). (6) 



Finally, given the approximation that the errors to [a/Fe] and [Fe/H] are uncorrelated, 
the error in the total metallicity is: 



0[M/H] = \/0"[Fe/H] 2 + 



0.638 * 10[ Q / Fe ] 
0.638 * 10[ Q / Fc l + 0.362 



07, 



t /FeY 



(7) 
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3. Fitting for the RGBB 



We jointly fit for the luminosity function of RG+RGBB stars in and near the RGBB 
using the following parametrization: 



N(I) = A exp 



B{I 

— I RGBB) 



N 



RGBB 



lit (T RGBB 



exp 



(I-I 



RGBB) 



2a 2 

^ RGBB 



(8) 



factoring out the normalization parameter A yields: 



N(I) = A< exp 



B(I - I 



RGBB) 



EWrgBB 



exp 



(I-I 



RGBB) 



^Crgbb 



(9) 



where A and B jointly define an exponential luminosity function for the underlying RG 
branch, the equivalent width EWrgbb is the ratio of the number of RGBB stars N RGBB to 
the number density of RG stars at the brightness of the RGBB, Irgbb is the mean brightness 
of the RGBB, and ctrgbb is the brightness dispersion of the RGBB. This methodology has 
previously been used elsewhere (Nataf et al. 2011afb ), but we will provide a stand-alone 
justification here. 

We fit in I because it is a more stable bandpass to blending and differential reddening 
than V, and because the absolute value of the derivative of the bolometric correction in 
I is much smaller than that in V for stars moving up to the RG branch, facilitating a 
more accurate comparison to models. We use I' = V— 0.6(5 — V^)— 0.4 for the WFPC2 
data, a relation obtained from the empirical calibration of T e g-logg-[Fe/H] ( An et aL||2007 ), 
calculated at the position of the RGBB. Vrgbb is obtained for both datasets by measuring 
the color of the RG branch at the position of the RGBB. CMDs, magnitude histograms and 
their corresponding best- fit N(I) probability density functions are shown for 47 Tuc, NGC 
362, NGC 1261 and NGC 7078 in Figures § [3l III U 



Measurement of the best-fit values for the parameters and their associated errors are 
done using a maximum-likelihood analysis to explore the parameter space via Markov Chain 
Monte Carlo (MCMC) - we do not bin the data before fitting the parameters. For each 
value of the parameters tested by the MCMC, we compute the log-likelihood £: 



N 



#(Ii/A,B,EWHGBB , orgbb, Irgbb) 



-N, 



obs 



(10) 



where iV bs is the total number of stars included in the fit. In each element of the chain, 
the parameter A is selected such that the integral of the function N(I) over the magnitude 
range is equal to N i, s . In other words, it is determined analytically for each combination of 
the other parameters rather than floated as a free parameter. 
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We can thus use the statistical identities valid for distributions marginalized to a single 
parameter, that parameter values with I > (£ max — 1/2) ar e within la of the best-fit, those 
with £ > i m ax — 2 are within 2a of the best-fit, and so on. This yields the reported 1-a 
statistical measurement errors for the parameters. A correlation diagram and the posterior 
distribution of parameters for 47 Tuc are shown in Figure [6j 

Since the number of stars expected in an interval (I,I+dI) is equal to N(I)dI, the 

statistical noise is Poissonian, thereby relating I to y 2 in the limit of a large number of 
datapoints: 

X 2 = -2i (11) 
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Fig. 2. — LEFT: CMD of 47 Tuc in ACS data, with the color-magnitude selection contours 
shown for the 2415 RG+RGBB stars and 545 RHB stars. The color of the RG branch at the 
RGBB is 1.03. The mean brightness of the RHB stars is Irhb = 13.09. RIGHT: Magnitude 
distribution of the RG+RGBB stars, Irgbb = 13.48 ± 0.01, EW RGBB = 0.32 ± 0.04. 
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Fig. 3.— LEFT: CMD of NGC 362 in ACS data, with the color-magnitude selection contours 
for 1059 RG+RGBB stars and 296 RHB stars shown. The color of the RG branch at the 
RGBB is 0.985. The mean brightness of the RHB stars is I RHB = 14.67. RIGHT: Magnitude 
distribution of the RG+RGBB stars, Irgbb = 14.41 ± 0.01, EW RGBB = 0.31 ± 0.06. 
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Fig. 4.— LEFT: CMD of NGC 1261 in ACS data, with the color-magnitude selection 
contours for 808 RG+RGBB stars. The color of the RG branch at the RGBB is 0.972. 
RIGHT: Magnitude distribution of the RG+RGBB stars, I RGB b = 15.63±0.01, EWrgbb = 
0.21 ± 0.06. The HB selection for this cluster can be found in Figure [8} 
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Fig. 5.— LEFT: CMD of NGC 7078 in ACS data, with the color-magnitude selection 
contours for 1403 RG+RGBB stars shown. The color of the RG branch at the RGBB is 
0.972. RIGHT: Magnitude distribution of the RG+RGBB stars, Irgbb = 14.24 ± 0.02, 
EWrgbb = 0.17 ±0.06. 

We thus kept, as part of our "gold" sample, all GCs that had a best-fit value of Nrgbb > 
10 and that were not known to be affected by severe patchy differential reddening. This 
sample on its own would be at risk of overestimating the expected value for EWrgbb, since 
it has a minimum value of the normalization. We construct a "silver" sample as follows. On 
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all the remaining clusters, we fit for the RGBB subject to the constraints: 

B = 0.72, (12) 

0.03 < cy rgbb < 0.09, (13) 
\Vrgbb ~~ Vrgbb ,predicted I < 0.2, (14) 
where the third constraint has a VRGss.predicted derived from the metallicity of the cluster 



and the V^-band apparent distance modulus to the cluster (Harris 1996, 2010 edition), and 
the best-fit relation for those two parameters obtained by a linear fit to the gold sample. 
The first two priors are absolute - no parameter space is explored outside the specified 
range. The third prior is slightly relaxed, other values of Vrgbb are explored, but with steep 
A£ penalties outside the specified range. To be included in the silver sample, GC CMDs 
had to first met one of three conditions, A^ RGj BB >measur ed > 5.0, N RG bb, predicted > 5.0, or both 
N rgbb, measured > 3.0 as well as A^GBB.predicted > 3.0. iV/jGBB.predicted is calculated from the fit 
to EWrgbb in the gold sample. We then also required that the error in the peak brightness 
of the RGBB be less than 0.05 mag. Larger values occur when the MCMC jumps between 
overdensities in the magnitude distribution - when it is not statistically clear which feature 
is the bump. We relaxed the a rgbb constraint to a maximum value of 0.12 mag rather 
than 0.09 mag for NGC 6316 and NGC 6440 due to their moderate differential reddening. 
GCs with severe patchy differential reddening such as NGC 6266 are not included. As the 
WFPC2 dataset has a heavier bias toward Disk/Bulge clusters, severe differential reddening 
proved to be a limiting criteria for a number of clusters. In total, there are 48 clusters in the 
gold sample and 24 clusters in the silver sample. 37 of the gold sample and 18 of the silver 
sample come from the ACS dataset, and the remainder come from the WFPC2 dataset. 



While we have used the parametrization discussed in this section before (Nataf et al. 
we recognize that it is a break from the great majority of the literature pertaining 
to the RGBB. Our investigation demonstrates that the RG luminosity function, the bright- 
ness peak, and normalization of the RGBB can be degenerate parameters (i.e. the errors of 
different parameters are correlated), and thus must be fit for concurrently rather than se- 
quentially. In light of this significant development in methodology, we independently discuss 
our parametrization for both the brightness and normalization of the RGBB. 



- 15 - 



NGC 104 (47 Tuc) Parameters 



420 

400 

< 380 

360 

340 
0.75 

m 07 

0.65 

0.6 
0.07 

m 0.06 

| 0.05 

b 0.04 

0.03 

0.5 

m 

g 0.4 

£0.3 

0.2 
13.5 
a, 13.49 
g 13.48 
~ 13.47 



Effl- 



13.46 =r 



Tl I | I I I I | I I I I | I I I I I I I I I |T 




cor=£k021 - 

■ .""'iisifi' 














r v 


cor 


=0 








cor= 




Ptt+f 



corF0.0S i 




/ i- 



— i — i — i — i — 

13.476 ±0.005 



-| — i — i — i — I — i — r 
■'COF<=- 0.543 




I I I I M I I I I I I I I 



cor=0.346 




cor,=Q.1.21 — 



334 




r I i i i i i "] I I I I I I I I i : 



i i i | i i _ 



0.33 ±0.04 



- 1 — i — I — i — i — i — I — i — r 
crc='-0.165 







±0.005 0.6 0.7 

B 



0.04 0.06 



RGBB 



0.2 0.4 
EW 

^ uu RGBB 



13.5 



RGBB 



Fig. 6. — Correlation diagram for a representative sample of points showing the distribution 
of the values of the RG+RGBB parameter limits from the MCMC for 47 Tuc. Parameter 
values within 4a of the best-fit value are shown in magenta, those within 3a in blue, those 
within 2a in green, and yellow for trials within la of the best-fit. Best-fit parameter values 
and the la errors are shown as a legend in the histograms, with correlations between different 
parameters shown as a legend in the scatter plots. Whereas the value of Irgbb is largely 
independent of the other parameters, that is not the case for EWrgbb- 
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3.1. The Brightness Peak of the RGBB 



Fusi Pecci et al. (1990) presented a method of measuring and interpreting the RGBB 



brightness that has since been broadly used. Their method was to log-integrate the lumi- 
nosity functions of cluster RG stars from both sides of the RGBB, and to measure the point 
at which the cumulative distribution function breaks with that of the two linear fits. The 
brightness is then compared to that of the zero-age horizontal branch (ZAHB) at the position 
of the RR Lyrae instability strip in V-band - the &Vrqbb = (Yrgbb — Vzahb) parameter. 
For this work, as well as our other recent works (Nataf et al. 2011a[b ), we use a maximum- 
likelihood method to fit for the parameters that has the advantage of being independent of 
bin size and of fitting all the parameters concurrently rather than sequentially. 

There are issues with the use of RR Lyrae stars as an anchor. First, ZAHB are much 
less likely to lie on the RR Lyrae instability strip in both metal-poor and metal-rich stellar 
systems. Additionally, in composite stellar populations such as dwarf galaxies, the Galactic 
bulge or indeed many massive GCs, the stars in the RR Lyrae instability strip may be 
biased toward a different subset of stars than the stars populating the RGBB. For example, 
the Galactic bulge RR Lyrae stars have a metallicity peak near [Fe/H]« —1.0 (Kunder & 



Chaboyer 2008 Pietrukowicz et al. 2011) compared to [Fe/H]« —0.1 for RG stars (Zoccali 



et al. 2008), and as such it would be unphysical to compare their luminosity to that of the 



stars in the RGBB. 

We suggest the use of two different comparative anchors for the brightness peak of the 



RGBB. First, following Cassisi et al. (2011), we compute wherever possible the difference in 
brightness between the RGBB and the main sequence turn-off, AV^f^ = (Vmsto—Vrgbb) ■ 
This parameter is more theoretically robust than AV^J^f , as it does not require assump- 
tions concerning the theory of post-main-sequence stellar evolution, such as neutrino energy 
loss. Moreover, for composite stellar populations, both the MSTO and the RGBB should 
be representative of the dominant population. We also compute the brightness parameter 
^Irgbb = (Irgbb — Irhb), where the mean brightness of the RHB is used as a benchmark, 
/-band is preferred due to the reduced evolutionary effects for the RHB in that bandpass 



(Girardi & Salaris 2001 Pietrzyhski et al. 2010). For heavily reddened systems such as 



some positions of the Galactic bulge (Nataf et al. 2011a), the RGBB may be measurable 



even when the MSTO falls at or below the photometric detection limit, whereas AI^q BB 
compares two populations at similar locations on the CMD. Additionally, the RHB is not 
only rigorously and thoroughly investigated in theory ( Girardi fc Salaris]|2001 ), it is the only 



Hipparcos-calibrated standard candle (Stanek & Garnavich 1998 Groenewegen 2008). For 



metal-rich populations, the RHB (as well as the MSTO) should share the RGBB's property 
of being representative of the numerically dominant population. 
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3.2. The Normalization of the RGBB 

Many prior studies of the RGBB in the literature have attempted to investigate star 
counts on the RGBB using the Reump parameter, which is the ratio of stars in the RGBB 
region Vrgbb — 0.4 < V < Vrgbb + 0.4 mag to that of stars with a brightness Vrgbb + 0.5 < 
V < Vrgbb + 1-5 ( |Bono et al.|2001[ ). We will not be using this parametrization, primarily due 



to its lower signal to noise ratio. Both characteristics can be traced to its arbitrary integration 
limits. The Reump parameter is also quite sensitive to photometric incompleteness. 

For typical values of EWrgbb and B, 0.3 mag and 0.72 mag -1 , the number of stars in 
the numerator of the RBump parameter, NReump, will be proportional to: 



NRsump = 0.3 + / exp 



o.i 



0.72(V - Vrgbb) 



dV 



-0.4 

NR B um P = 0.3 + 0.81 = 1.11. (15) 
Meanwhile for the denominator: 

'1.5 



DRsump = / exp 

'0.5 



0.72(V- V RG bb) 



dV = 2.10. (16) 



Only (0.3/1.11) ~27% of the stars in the numerator come from the RGBB stars - the rest 
measures the underlying RG branch. To see the consequences of this, consider the example 
of a well-populated CMD that is expected to have ~100 RGBB stars, more than the number 
in 66 of our 72 clusters. Such a system would then also be expected to have 270 stars RG 
in its numerator, and 700 stars in its denominator. The signal to noise ratio of the RGBB 
would then only be a meager 100/sqrt(100+270+700) ~ 3.0. 

As can be seen from the calculation above, a GC without a single RGBB star but with 



a typical RG luminosity function would have an Rsump value of 0.81/2. 10~0. 39. Bono et 



al. (2001) list RBump for 47 Tuc as having a value of (0.63 ± 0.05) - a 4.8<r detection. By 
contrast, using our parametrization, we measure N R q BB = (122.3 ± 14.2) for 47 Tuc - an 
8.6cr detection. 

The situation will then be much worse for the vast majority of GCs that have less well- 
sampled CMDs, and/or for the metal-poor clusters that have a lower value of EWrgbb- 
Further, a signal of <3.0<r is just for the zeroth order existence of the RGBB. The signal 
will be much lower if one investigates first-order effects such as gradients due to age, he- 
lium and metallicity. Fundamentally, RBump is a composite parameter of the parameters 
B and EWrgbb, with a heavy bias toward B. We consider both these parameters to be 
independently interesting, and argue they should be fit as distinct parameters. 
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The two normalization parameters we introduced in Nataf et al. (2011a) and study 
in further detail here mitigate this issue. The EWrgbb and f^Q BB , are the number of 
RGBB stars respectively normalized by the number density of RG stars per magnitude at 
the brightness peak of the RGBB, and the number of HB stars. As the RGBB is only 
observed as an excess over the continuum luminosity function of the RG branch, its best-fit 
normalization will always be degenerate with the parameters A and B. We reduce the impact 
of this degeneracy by including as many stars in the fit as we can while simultaneously leaving 
out contamination from the SGB, the HB, the AGB, and foreground disk contamination 
where present. We do not discard the statistically meaningful stars that are either brighter 
than the RGBB by more than 0.4 mag, between 0.4 and 0.5 mag fainter than the RGBB, or 
fainter than Irgbb + 1-5 but still markedly brighter than the SGB. 

Broad applicability ought to be a high priority for the definition of any astrophysical 
parameter, and RBump does not generalize well to composite stellar systems. In our investi- 



gation of the Galactic bulge RGBB (Nataf et al. 2011a), we did not integrate up to stars 1.5 



mag fainter than the RGBB. That region of the CMD is heavily mixed with foreground disk 
stars and bulge SGB stars. The integration limits need to be flexible in order to account for 
the diversity of stellar populations in which the RGBB is observable and will be observable 
in the future. 

Lastly, photometric incompleteness will sometimes be a concern. It is true that both of 
the catalogs used in this work have artificial star tests confirming completeness on the RG 
branch, but there is value in having a methodology that could generalize well to other kinds 
of catalogs. For any smooth photometric completeness function, the parameter EWrgbb 
has the advantage of not incorporating photometric incompleteness as a systematic error, 
since RGBB stars will have the same detection probability as the RG stars at the same 
brightness. The systematic effect will be limited to reducing the value of B, an unfortunate 
but contained issue. The case is different for the RBump parameter. Since the fainter stars 
in the denominator will be less frequently detected, photometric incompleteness will end up 
masquerading as a stronger normalization for the RGBB. 



3.3. The Continuum Distribution of RG Stars 

Fitting an exponential law to the number distribution of RG stars as a function of 
magnitude has the simple physical explanation that it corresponds to a power-law as a 
function of luminosity. An exponential fit to the number counts also corresponds exactly 
to a linear fit to the log of the number counts. Previously, the RG continuum distribution 
has been modelled by fitting a linear relationship between the log of the cumulative number 



-19- 



counts and magnitude (Fusi Pecci et al. 1990 Zoccali & Piotto 2000). We prefer to use 
the number density functions rather than the cumulative density functions because it is 
straightforward to calculate the errors for the former. Cumulative distribution functions do 
not have straightforward error calculations because adjacent bins have correlated number 
counts. 



4. The Main-Sequence Turnoff 

We fit for the MSTO in all of the ACS GCs for which we have a measurement of the 



RGBB. We do not use the MSTO measurements of Marin- Franch et al. (2009) since those 
are reported in the (F606,F606W — F814W) absolute magnitude plane rather than in the 
(I,V — I) apparent magnitude plane required by our investigation. 

To measure the MSTO of each cluster, we first fit 2nd and 3rd degree polynomials to 
the upper main-sequence of each cluster, with color being a function of magnitude. The 
fits typically cover a luminosity range between 0.3 mag brighter than the MSTO and 0.7 
mag fainter, the boundaries are selected to comfortably include the MSTO but to exclude 
regions of the CMD that would distort the fit, such as the subgiant branch. 3-cr outliers are 
recursively removed from fits with the fits then recomputed, though we removed 2-cr outliers 
in NGC 6171 and NGC 6624 due to their thicker main-sequences. For both polynomials, we 
take the bluest point on the best-fit curve as an estimate for the MSTO, and we report the 
average of the two values as our measurement. 

The values obtained by the 3rd and 4th order polynomials had an average discrepancy 
of 0.01 mag in V and 0.0002 mag in (V — I). Unfortunately, it is difficult to quantify 
the errors in these fits due to the existence of a few systematics. For example, we expect 
some contamination from binary stars, though the most egregious examples of those are 
left out of the fit by our removing of outliers. We show the results of our method for two 
representative GCs, NGC 104 (47 Tuc) and NGC 1261 in Figure (7) The summary of our 
MSTO measurements in Table [U 
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Fig. 7.— TOP: CMD of the relatively metal-rich, disk/bulge cluster NGC 104 (47 Tuc) 
shown, with the measured position of the MSTO, (V — I,V) = (0.696,17.710) delineated 
by the white point. Right panel shows the CMD zoomed in, with our best-fit polynomial 
to the main sequence shown as the white line and the measured position of the MSTO as 
the black point. BOTTOM: Same as in the top panel, but with the CMD shown for the 
metal-intermediate, outer halo cluster NGC 1261. The position of the MSTO is (V — I, V) = 
(0.584,20.115) 



5. The Horizontal Branch Stars 

We estimate the number of HB stars in all GCs for which we have a measurement of the 
RGBB. We also compute the mean brightness of the RHB for 31 of the clusters for which the 
RHB is observed to be a well-populated, visually distinct component of the CMD. The error 
in the mean brightness, o~i RHB , is taken to be the standard error in the mean, oj \/{N RHB ). 
Two examples of RHB selection cuts, NGC 1261 and NGC 7089 (M2), are shown in Figures 
[8] and [9j Additionally, the examples of 47 Tuc and NGC 362 are shown in Figures [2] and 
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|3j The RHB is conservatively selected by drawing a box where a clear RC is visible. GCs 
without clear RCs do not contribute to the AI B q BB statistic. 

Due to the fact that the HB populates a specific region of the CMD, it is relatively 
straightforward to count up the number of HB stars. However, as is clearly discernible 
in Figures [8] and |9j blue stragglers, background (or foreground) contamination stars, and 
AGB stars are sometimes photometrically indistinguishable from HB stars. Fortunately, 
the intersection of those populations with that of the HB on the CMD never totals more 
than a few percent of the HB population. Since the uncertainty in the number of RGBB 
stars is typically ~10%, the small systematic uncertainty in the HB number counts does not 
contribute to the error budget of the parameter f^Q BB . 



Zoccali & Piotto (2000) counted the HB stars in 26 of the GCs observed in the WFPC2 



survey (Piotto et al. 2002). Our investigations have 8 CMDs in common. We measured 168 
HB stars in NGC 1904 to their 177, 145 HB stars in NGC 5634 to their 146, 529 HB stars 
in NGC 5824 to their 520, 302 HB stars in NGC 6139 to their 299, 34 stars in NGC 6235 to 
their 35, 133 HB stars in NGC 6284 to their 133, and finally 365 HB stars in NGC 6356 to 
their 370. The level of disagreement is thus of order 2%, miniscule compared to the typical, 
~10% error in Nrgbb, or even the typical ~7% Poisson error in Nhb- 

We recognize that there may be a small bias in our measurement of J^gbb due t° the 
fact these HST data are taken toward the cores of GCs. Due to chemically-distinct multiple 
generations (D'Ercole et al. 2008) and dynamical relaxation (Leigh et al. 2011), different 
phases of stellar evolution should have slightly different occupation ratios at differing core 
radii. It is difficult to assess the impact of these effects at this time since they are both 
rapidly evolving fields. 



The HB characterization for the clusters NGC 6388 and NGC 6441 are somewhat mod- 
ified due to their distinct morphologies. These are discussed in Section 9.3 
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Fig. 8.— ACS CMD for the intermediate-metallicity GC NGC 1261. The 231 HB stars are 
shown by open circles. The mean brightness for the RHB stars is Irhb = 15.94. Some MS 
stars have been removed from the figure to reduce image size. 
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14 - 



16 



18 



20 - 




-0.5 







0.5 



V - I 



Fig. 9.— ACS CMD for the low-metallicity GC NGC 7089 (M2). The 720 HB stars are 
shown by open circles. Some MS stars have been removed from the figure to reduce image 
size. 



6. Results: The Brightness and Color of the RGBB 

We have measured the brightness of the RGBB in 55 of the clusters from the ACS 
dataset and 17 from the WFPC2 dataset, for a total of 72 measurements. All of the ACS 
clusters have measured values of the MSTO positions and thus values of AV^sto an d 
A(V - I)ffi S TO- The relative brightness between the RGBB and the RHB, AI%g£ B , for 22 
of the ACS clusters and 9 of the WFPC2 clusters, for a total of 31 measurements. All the 
measurements discussed and used in this section are listed in Tables 4 and 5. The GCs 
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NGC 2808, NGC 5286, NGC 6388 and NGC 6441 are not included in the fits due to their 
anomalous RGBB properties. 



The dominant empirical trend, previously observed by several investigations (Fusi Pecci 



et al. 


1990 


Zoccali et al. |1999 I 


Hello et al. 2003 


Monelli et al. 


2010 


Nataf et al. 


2011a 


), is the c 



increasing metallicity. The Mv,rgbb increases by ~1.6 mag as the metallicity increases 
from [M/H]w -2.1 to [M/H]« -0.1. The change in brightness relative to the MSTO is 
shallower, ~1.1 mag over the metallicity range, due to the fact the MSTO also gets fainter 
with increasing metallicity. For AI^q BB = (Irgbb — Irhb), we find a variation of 1.1 mag 
over a metallicity interval of ~1.3 dex. 

We compute linear fits for all three of the brightness variables. For the brightness 
Mv,rgbb, a fit weighted by the statistical error in the brightness measurements has a 
X 2 =3785. The clear interpretation is that the dispersion due to errors in the input metallic- 
ities, apparent distance modulus in V, and undiagnosed second parameters are substantially 
larger than the statistical error in the measurement of the brightness. We adjust the errors 
using the following prescription: 



crM. 



v,rgbb 



2' 



oM- 



V,RGBB 



+ 



V.RGBB 



^[M/H] 2 + 5aV R GBB 2 , 



(17) 



d[M/H] 

where 8<jMv,rgbb is the noise added due to undiagnosed second parameters. The fit obtained 
is My rgbb — (0.600±0.013) + (0.737±0.024)([M/H] + 1.110). A value 6aM VtRGBB = 0.077 
mag is needed to yield a fit with \ 2 — 66 (68 measurements and 2 parameters). This scatter 
could also be due to errors in the values of the V^-band apparent distance modulus summa- 



rized by Harris (1996, 2010 edition), as well as the fact that many of the distance estimates 
come from different methods, rendering the list of distance moduli used heterogeneous. The 
scatter in Mv,rgbb could also be due to an additional scatter of ~ (0.077/0.737 = 0.104) 
dex -1 in the metallicity [M/H] above that which is assumed in this work. Another possibility 



are variations in age or initial helium abundance (Cassisi & Salaris 1997). 



For the brightness relative to the MSTO, we obtain AV$g£g = (3.565±0.012) + 
(-0.549±0.023)([M/H] + 1.152). The intrinsic scatter required to have \ 2 =49 is 8AV$gfg = 
0.072 mag. Age could be the source of this extra scatter. Stellar evolution models predict 



that older clusters should have larger values of AV^f-f^ , at a rate of ~0.05 mag/Gyr (Cas- 



sisi et al. 2011). Both the RGBB and MSTO become fainter with increased age, but the 
MSTO becomes fainter faster. The extra scatter in AV^sto can thus be entirely explained 
by an age scatter of ~1.5 Gyr for GCs with respect to the age-metallicity relation of GCs. 
It is noteworthy that 6AV$gfg =0.072 mag is almost as large as 5 a My rgbb = 0.077 
mag. This suggests that either the apparent distance modulus is accurate to a level of 
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v/0.077 2 - 0.072 2 « 0.03 mag, or that the errors in the V^-band apparent distance modulus 
correlate with the sources of variation in the position of M^msto- The difference in (V — I) 
color decreases with increasing metallicity. As metallicity increases, both the MSTO and 
RGBB become redder, but the MSTO becomes redder faster. 
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Fig. 10. — In all 4 panels, blue points are the gold sample measurements in the ACS clusters, 
red points are the gold sample measurements in the WFPC2 clusters, and green points come 
from the combined silver sample. TOP-LEFT: ^V^qbb = (Ymsto ~ Vrgbb) for all the 
ACS GCs. TOP-RIGHT: A(V - I)%Sbb = (V - I)msto - (V - I) RGB b for the ACS 
GCs. BOTTOM-LEFT: M-v,rgbb for all clusters with an RGBB measurement, using the 
V^-band apparent distance modulus from Harris (1996, 2010 edition). BOTTOM-RIGHT: 
^Irgbb = {Irgbb — Irhb) from 31 GCs for which we measured the RHB mean brightness. 
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For AIfgl§ B , we obtain AI$J§ B = (0.123±0.018) + (0.852±0.045)([M/H] + 0.640). The 
measured intrinsic scatter is 5AI B q BB = 0.051 mag. 5AI^q BB is smaller in quadrature than 
SAVfiQBB by a value of 0.051 mag, so it may appear to be a more stable variable. However, 
it is calculated from a sample of 28 rather than 51 GCs, and over a smaller metallicity 
range. Moreover, it has a hidden selection bias. Whereas all metal-rich clusters have a red 
component to their HB, only some of the intermediate-metallicity ([Fe/H]~ —1.5) GCs do. 



This is due to the effect of second parameters, possibly age (Lee et al. 1994 Dotter et al. 
201TJ. 



7. Results: The Number Density of RGBB Stars 



We measure fgg BB and EW RGBB in 55 of the ACS GCs and 17 of the WFPC2 GCs for a 



total of 72 measurements. The measurements are shown in Figure 11, and all measurements 
discussed in this section are listed in Table 6. The GCs NGC 2808, NGC 5286, NGC 6388 
and NGC 6441 are not included in the fits due to their anomalous RGBB properties. 

The dominant empirical trend is the increasing number counts of the RGBB with in- 
creased metallicity. There are two factors involved. The first is that the RGBB gets more 
prominent relative to the underlying RG branch at increased metallicity. As [M/H] is in- 
creased from —2.0 to 0.0, EWrgbb increases by a factor of ~2.7. Further, there is the 
additional effect that evolution is slower further down the RG branch, amplifying the first 
effect. As [M/H] is increased from —2.0 to 0.0, f B c BB increases by a factor of ~8.0. 

Some care must be taken in computing a linear fit for these parameters. There are 
significant statistical errors due to the error in the measurement ofNucBB, as well as expected 
fluctuations due to hidden second parameters such as variations in age, CNO abundances, 
initial helium abundance, and other factors, since the strength of the RGBB will not be a 



function of metallicity alone (Nataf et al. 2011a). For example, a fit for f B Q BB with respect 



to metallicity weighted purely by the statistical error measurements yields \ 2 = 96.8 for 68 
measurements and 2 parameters. It is clear that whereas metallicity is the first parameter 
of RGBB strength, it is not the only parameter. 

We add, in quadrature, a systematic noise to the error Scrf B Q BB with the following 
prescription: 

tHB 12 



_/ fHB 2 _ _ pHB 2 , 
a J RGBB — J RGBB <~ 



RGBB 



C[M/H] +^crf RGBB , (U 



_d[M/H]_ 

And we compute a weighted least squares using the combined error o' f B Q BB , 

Wi = l/(a7g? M ) 2 , (19) 
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where we adjust the value of Saf B Q BB until we obtain a fit with x 2 = 66. The analo- 
gous procedure is performed for EWrgbb- We thus measure f B Q BB = (0.111±0.005) + 
(0.109±0.011)*([M/H] + 1.273), with a second parameters noise value of 6cf%g BB = 0.018. 
Similarly we fitted EWrgbb = (0.248±0.010) + (0.121±0.018)*([M/H] + 1.134). X 2 = 62.1 
for 68 measurements and 2 parameters, implying a perhaps surprising lack of evidence for 
hidden second parameters. Nataf et al. (2011a) found a predicted variation of cIEWrgbb/ 'dt 
= 0.008 mag Gyr -1 for a 10 Gyr old population with [M/H]=0, an effect which could likely be 
too small to infer given the large statistical errors. More theoretical investigation is needed 
to ascertain whether the value of EWrgbb is slowly varying with age across the range of 
metallicities probed in this work. Moreover, our fits may have too great a degree of freedom. 
As we will show in the next section, the RG luminosity parameter B may in fact be a con- 
stant or a quasi- const ant of stellar evolution. If this is confirmed by subsequent theoretical 
investigation, then future fitting would need to take that into account, and the parameter 
space of values for EWrgbb would shrink accordingly. 
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Fig. 11. — For both panels, blue points are measurements in the gold sample of ACS clusters, 
and red points are measurements in the gold sample WFPC2 clusters. The points from the 
silver sample are shown in green. The best-fit line is shown for both relations. TOP: 
EWrgbb for 63 GCs as a function of metallicity. BOTTOM: The fraction of RGBB stars 
to HB stars clS db function of metallicity. 



8. Results: Other Parameters 

We briefly discuss the other measured stellar evolution parameters, the exponential slope 
of the RG luminosity function B, and the magnitude dispersion of the RGBB ctrgbb- The 
parameters discussed in this section are listed in Table 7. We do not include GCs from the 
silver sample in our fits for these two parameters, as these GCs had these two parameters 
constrained to match the distribution in the gold sample. Therefore, the fits in this section 
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are done purely on the 44 measurements in the gold sample, with NGC 2808, 5286, 6388 
and 6441 removed as before. 
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Fig. 12. — TOP: B, the exponential slope of the RG luminosity function, as a function of 
[M/H]. BOTTOM: The measured magnitude dispersion of the RGBB. 

We find that B has no significant dependence on metallicity. For a least squares weighted 
by the errors in the measurement of B and using only the measurements in the gold sample, 
we obtain 5=(0.715±0.006) + (0.008±0.012)([M/H] + 1.105). x 2 =47.2 for 44 measurements 
and 2 parameters. The slope is detected at the ~0.6<r level - it is not significant. A fit to 
the weighted mean value B = 0.719 yields a \ 2 =47.7. In light of the potential systematics 
present such as varying amounts of disk contamination in clusters, we argue that there is no 
convincing evidence for a relation with metallicity. The mean of the measurements weighted 
by the errors in B is (0.715±0.006). 
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This behavior in B is predicted by stellar models, as a straightforward consequence of 
the relation between the total luminosity of a star on the RG phase and the mass of the 



He-core (Paczynski 1984; Castellani et al. 1989). The prediction of Castellani et al. (1989), 
that B = 0.74 ± 0.04 across the age and metallicity range spanned by the Milky Way GC 
system is confirmed by our investigation. 

We also measure the relation for the magnitude dispersion of the RGBB. (Trgbb — 
(0.051±0.003) + (0.017±0.006)*([M/H] + 0.709), with X 2 = 32.8 for 44 measurements and 
2 parameters. Stellar evolution models predict that more metal-rich stellar systems should 



have broader bumps (Fusi Pecci et al. 1990 Cassisi & Salaris 1997 Di Cecco et al. 2010). 



However, broader RGBBs can also be obtained from differential reddening, blending, and 
from multiple populations. Due to the many contributing factors, we think that while 
crgbb may be useful in interpreting the bumps of specific GCs, one should be cautious in 
interpreting the global relation. 



9. Interesting Clusters: NGC 2808, 5286, 6388, and 6441 

We comment on the interesting anomalies we measure in the RGBBs of the GCs NGC 
2808, 5286, 6388 and 6441. 



9.1. NGC 2808 



NGC 2808 is known to have at least three main sequences, from which |Piotto et al. 
(2007) estimates two helium-enhanced populations, each with ~15% of the cluster stars. 



Their inferred enhancements are AY~0.05 and 0.12. Due to theoretical expectations of the 
effects of helium on the RGBB (Cassisi & Salaris 1997 Nataf et al. 2011a[b ), one should 
expect the RGBB to be slightly brighter and less populated, though the effect will not be very 
strong if 70% of the stars are first-generation, and the fact second-generation stars are also 
enhanced in total (C+N+O) ( Salaris et al.||2006 ). The bigger impact will be on the shape of 
the RGBB. We test for skew-normal distributions on all 5 of the ACS clusters that have at 
least 50 RGBB stars and that do not have a RHB mixed with their RG+RGBB branch: 47 
Tuc, NGC 1851, NGC 2808, NGC 5927, and NGC 6624. The A X 2 values are 1.50, 0.08, 3.49, 
0.66, and 0.53 respectively. Only NGC 2808 exhibits a strong detection of a skewed RGBB. 
Its parameters change from (Vrgbb, <?rgbb, EWrgbb) = (16.235, 0.092, 0.347) to (16.219, 
0.112, 0.303) with a strongly negative skew of — 0.924^q q33. The negative skew is exactly 
what one would expect if there were a relatively small number of brighter (helium-enhanced) 
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RGBB stars. Moreover, it also contrasts to the expectation from models ( Cassisi et al.|[2002 



Nataf et al.||2011b[ ) that the RGBB of a single-metallicity, single-age population be positively 



skewed, i.e. with its mode at its bright end and a long tail to fainter luminosities (and thus 



higher values of magnitude). We show the CMD and magnitude histogram in Figure 13 




Fig. 13.— LEFT: CMD of NGC 2808 from the ACS sample, zoomed in on the RG branch 
at the location of the RGBB. RIGHT: Magnitude histogram of 3308 RG+RGBB stars with 
fit. Unlike other well-populated clusters, the RGBB is much better fit by a skew-normal 
distribution than by a standard normal distribution, consistent with findings that the cluster 
has an extreme helium-enhancement subpopulation (Piotto et al. 2007). 



Consistent with the measurement of the third moment (skewness) of the RGBB's mag- 
nitude distribution is that of the second moment (dispersion). At the cluster's metallicity 
[M/H]= -0.94, the predicted width from the total GC sample is a RGBB = 0.047±0.003. 
We measure urgbb = 0.092 ± 0.013 without a skew and 0.112±0.015 with a skew. The 
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estimated differential reddening in the cluster is AE(B-V) ~0.02 mag (Piotto et al. 2007), 



corresponding to AAj ~0.036 mag, assuming a standard extinction law ( Cardelli et al.|19 89). 
This leaves a significant source of dispersion of 0.07-0.10 mag due to multiple populations. 



Bragaglia et al. (2010) predicted a broadening of the RGBB in clusters with He-enhanced 



populations with AY > 0.10, due to the magnitude separation between the RGBBs of the 
first and second generations. 

The second and third moments of the NGC 2808 RGBB are in a rather startling agree- 
ment with stellar evolution predictions and observations of other stellar populations of the 
cluster. A clear prediction is that the brighter RGBB stars should be oxygen-poor and 
sodium-rich relative to the fainter RGBB stars, since those are thought to be the most 
helium-enriched (Piotto et al. 2007; Valcarce & Catelan 2011). With a population of ~170 



RGBB stars in this CMD alone, NGC 2808 may be one of the few Galactic GCs in which 



this experiment is feasible. Bragaglia et al. (2010) did detect a global brightness difference 
of ~ 0.044 ± 0.042 mag difference in V for 1368 stars from 14 GCs characterized as either 
primordial, or intermediate/extreme based on their [Na/Fe] abundances. 



9.2. NGC 5286 

NGC 5286 may be displaying a split RGBB. We show the CMD and magnitude his- 



togram in Figure [14} There are two peaks, one at Vrgbb = 16.287 and a second, smaller 
peak ~0.2 mag fainter. On its own, we do not consider the second peak to be individually 
compelling, but we do consider it worthy of mention in light of other issues. The value 
of MyjiQBB for the brighter, more populated peak is ~0.15 mag brighter than the expec- 
tation from the fit to all the GCs, and AV^q^q = 3.87, is larger than the expected 3.72 
by a similar amount. The number counts are also low. At [M/H]= —1.43, the predicted 
value of EWrgbb from the fit is 0.212±0.011, whereas the measured value in this GC is 
0.123±0.0136. Nearly one half the RGBB stars are "missing". With a sample of 23.4±6.6 
RGBB stars, a Poisson fluctuation of the required amplitude is very unlikely. 

The right panel of Figure [14] presents a simple solution - one third the RGBB are located 
in the second overdensity ~0.2 mag fainter than the RGBB. This secondary overdensity is 
also seen in the more comprehensive (B — V,V) CMD of NGC 5286 from ground-based data 



shown in Figure 8 of Zorotovic et al. (2009), where it is quantified as a ~4a effect. The 
weighted mean of the brightnesses would yield the approximate expected values of AVb G S bb 
and Vrgbb, and the sum of their normalization would do likewise for EWrgbb and /rgrb- 

One solution is for this cluster to be an extreme member of the class formed by NGC 
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1851 and NGC 6656 ( Milone et al.|2011 ). Those two GCs do not show the behavior expected 
from a spread in helium, but spectra demonstrate variations in heavy elements such as iron 
and yttrium. The two peaks are matched by theory if two thirds of the stars (the brighter 
RGBB) are in the first generation, and the remainder are in the second, [Fe/H]-enhanced 
second generation. Zorotovic et al. (2010) inferred the [Fe/H] of the RRc stars from their 
Fourier components. The [Fe/H] values are -1.90, -1.90, -1.89, -1.89, -1.84, -1.77, 
— 1.69, —1.68, —1.67, —1.66, —1.61, and —1.07. The cumulative distribution is shown in 
Figure 15, This is quite consistent with two peaks, one near —1.85 and the other near 



— 1.65, with one outlier with [Fe/H]= —1.069. The mean of these measurements excluding 
the outlier is [Fe/H]= — 1.771±0.036. Only one of the stars is near this mean. 




Fig. 14.— LEFT: CMD of NGC 5286 from the ACS sample, zoomed in on the RG branch 
at the location of the RGBB. RIGHT: Magnitude histogram of RG+RGBB stars with fit. 
A secondary overdensity is seen at a magnitude ~0.2 mag fainter than the main RGBB. 
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We summarize the lines of evidence for this cluster being an extreme analog of NGC 
1851 and NGC 6656: 

1. The RGBB measured by the fit is significantly brighter (0.10-0.15 mag) than expected 
from the relation derived from the GCs studied in this work. 

2. The RGBB measured by the fit is significantly less numerous (50%) than expectated 
from the fit derived to the GCs studied in this work. 

3. Both of the first two concerns are resolved if one takes the RGBB to be the weighted 
sum of the peak found by the single-peak fitting code and the smaller peak that is 0.2 
mag fainter. 



4. the RRc stars studied by Zorotovic et al. (2010) show a double peak in their inferred 



[Fe/H] values, on opposite sides of and similarly displaced from the photometric metal- 
licity, with both corresponding well to the two putative RGBB peaks. 




Fig. 15. — The cumulative distribution of [Fe/H] for the RRc stars in NGC 5286 derived from 
Fourier coefficients by Zorotovic et al. (2010). The values appear clustered at [Fe/H]= —1.85 
and [Fe/H]= -1.65. 
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Unfortunately, there does not appear to have been much spectroscopic investigation of 
this GC, so we cannot know for certain at this time. We check the online version of [Harris 



| (1996, 2010 edition) and find that there is only one reference for the metallicity, and it is 
a photometric metallicity. A photometric metallicity would necessarily give the weighted 
mean of any metallicity spread, with the color width in the stellar populations being well-fit 
by differential reddening or photometric noise. We have found an older reference in the 



literature (Rutledge et al. 1997), which inferred metallicities for 52 Galactic GCs using Call 



triplet equivalent widths. Figure 7D shows some sccatter in the 17 EWs measured for RGs 



in NGC 5286, but further investigation would be required to confirm. Rutledge et al. (1997) 



reported [Fe/H]= — 1.70±0.03 as the most probable value on the metallicity scale of Zinn & 



West (1984) 



While none of the four lines of evidence listed above may be individually convincing, 
their union constitutes a compelling case to obtain high-dispersion spectra of RG stars in 
this cluster. 



9.3. NGC 6388 and NGC 6441 

We find a few peculiarities in the GCs NGC 6388 and NGC 6441. Their relative mea- 
surements are not consistent with their identical spectroscopic metallicities, the brightness 
of the RGBB indicates the distance may be underestimated, and their low number counts 
are consistent with the presence of an extreme, helium-enhanced population. 

We first state the different HB calibration selection used for these two GCs. For these 
GCs, a quantitatively significant portion of the red end of their RHB merges with the RG 



branch. Due to the fact these are two of the most interesting GCs in the Galaxy (Yoon et 



al. 2008), it is critical to adapt our method to get these right. We first sum the number 
of point sources toward the regions of the CMD that are clearly dominated by the HB. We 
then fit for a second Gaussian in the combined RG+RGBB+RHB branch to measure the 
red end of the RHB. The total HB population is then the sum of the HB stars counted in the 
rest of the CMD and the best-fit normalization value of the RHB component along the RG 
branch. We take the weighted mean (by number counts) for the brightness of the RHB. It is 
necessary to do a double-Gaussian fit if only to have a proper fit of the RG+RGBB+RHB 
branch. Without doing this the parameters of the RGBB would be severely compromised. 



The procedure is visually summarized in Figure 16 
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Fig. 16.— LEFT: CMD of NGC 6441 in ACS data, with the color-magnitude selection 
contours for 5776 RG+RGBB+RHB stars and 1280 pure RHB stars shown. The color of 
the RG branch at the RGBB is 0.985. The weighted mean brightness of the 1280 RHB 
stars in the pure RHB box and the 433 mixed with the RG+RGBB stars is Irhb = 16.38. 
RIGHT: Magnitude distribution of the RG+RGBB+RHB stars, I RGBB = 16.80 ± 0.01, 
EW RGBB = 0.22 ± 0.02. 



The brightness parameters for these two GCs are not compatible with their reported 
metallicities. [M/H] 64 4i-[M/H] 6388 = 0.01+0.08 dex flCarrettaet al.|2009b[ ), but AI^ B \ eui - 
AIrgbbW388 = 0.109+ 0.0 1 7, perhaps indicating an actual difference in metallicity of ~0.13 
dex. Both clusters have very faint values of Mv,rgbb even as they have values of AV^^g 
that agree with the global trends to within ~0.05 mag, implying an error in either the dis- 
tance or reddenings to the clusters. NGC 6388 has a Vrgbb value 0.32 mag fainter than 
that predicted by the fit, and the deviation is 0.41 mag for NGC 6441. These deviations are 
far larger than could be reasonably attributed to variations in age or to errors in metallicity, 
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so we conclude that a combination of the distance and reddening to the clusters listed in 



Harris (1996, 2010 edition) are underestimated. The apparent distance modulus is obtained 



from RR Lyrae measurements (Pritzl et al. 2001, 2002). If the RR Lyrae stars are helium 



enhanced, then they will be brighter than that predicted by the standard [Fe/H]-My relation 



for RR Lyraes (Caloi & D'Antona 2007), leading to a severe underestimate of the distance. 



We argue this to be the case here. 

Enhanced helium enrichment may play a role in the RGBB star counts for this cluster. 
The fncBB derived from the relation to all the GCs at [M/H]= -0.295 is 0.217±0.012, 
whereas we measure 0.156±0.017 and 0.151±0.016 for NGC 6388 and NGC 6441. A ~30% 



deficiency is detected at the level of 2.9<r and 3.3<r respectively. Caloi & D'Antona (2007) 



argued that at 15% of the stellar content of NGC 6441 had to be extremely helium enhanced 
(AY > 0.1) in order for stellar models of HB evolution to match observations of the clusters, 
in particular the well-populated RR Lyrae instability strip and blue horizontal branch (BHB) 
at the high-metallicity of the cluster. If that is the case, not only would the helium-enhanced 
stars have shorter RGBB lifetimes, but their characteristic magnitude would be at a different 
location on the RG branch. The low values obtained for RGBB star counts could thus 
be interpreted as being due to only most, rather than all, of the RGBBB stars having a 
brightness at or close to that of the measured peak. If ~30% were significantly brighter, and 
fell between the RGBB and the RHB or perhaps even in the RHB region of the CMD, they 
would not be captured by our measurement. A similar scenario may be at work for NGC 
6388. 

Both GCs have values of orgbb that are larger than expectations from the fit to all the 
GCs. However, it is difficult to interpret these excesses due to the large uncertainty in the 



differential reddening toward these GCs (Yoon et al. 2008). 



Unfortunately, disk contamination could in principle corrupt the measurements in NGC 



6441, as can be clearly seen in Figure [16} We experimented with various selections for the RG 
branch, shifting the limits for both color and magnitude, and our measured values of /robb 
did not change by more than 1%, nor was the brightness peak shifted. The measurements 
also remained the same when we removed the inner half of the GC stars, to test for effects 
due to photometric noise. As well as we can test with the available data, the measurement 
of a discrepancy appears robust. 
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10. Comparisons to the Milky Way Bulge and M32 



The first measurement of the RGBB in the Milky Way bulge was discussed in Nataf et 



al. (2011a). The detection was later confirmed by Clarkson et al. (2011) and Gonzalez et al. 



al. 


(2011a 


(2011b 


). S 



cRC 



, ulge = 0.201±0.012, 

A-fflGRB I Bulge =0.737±0.012, where we use "RC" for red clump to refer to the RHB of these 
populations for consistency with the literature. These values are a little different than those 



reported by Nataf et al. (2011a) due to improved quality criteria: we select cleaner sightlines. 



The parameters are measured toward the sightline (l,b) = (0, — 2), where the brightness 
dispersions are (Jrgbb = 0.220±0.010 and ctrgbb = 0.241±0.003. As discussed by |Nataf 



& Udalski (2011), of the available OGLE-III sightlines toward the "triaxial bar" compo- 
nent of the bulge, this is the sightline with the smallest quantity of observed geometric 
dispersion, and therefore the minimum degeneracy between the different components of the 
RC+RG+RGBB luminosity function. 



The measurement for the dwarf elliptical M32 comes from the investigation of Monachesi 



et al. (2011). Within their field, imaged by HST, they report a value of Nrc = 1422.8, and 
Nrgbb = 219 ± 51, for /J§bbIm32 = 0.151±0.036. The difference in brightness reported is 
M32 =0.56 ±0.13. Their Figure 12 shows that the two features are at equal or 
very nearly equal colors, so we adopt the value A/Jjq BB |m32 =0.56 ±0.13. 



These are the parameters we show in Figure [T] and that we assume for discussion through- 
out this work. In particular, we find that both systems have lower values of /J^Jbb than that 
expected from the Galactic GC system, and that the bulge has a lower value than expected 
for A/ BGBB | Bulge once composite metallicity effects are taken into account. Because these 
two measurements were obtained with different instrumentation, different methodology, and 
other different systematics such as crowding, their comparable deviation from the relations 
for Galactic GCs is independently derived, and may be due to a similar evolution. 



11. Application: Empirically-Motivated Prediction of the Galactic Bulge 

RGBB Properties 



The Galactic bulge is a complex stellar population for which different analyses lead to 



different results for parameters as fundamental as age (e.g. Bensby et al. 2011 Clarkson 



et al. 2011). In this section, we compute what the Galactic bulge RGBB should look like 



relative to the RC toward two distinct sightlines given the assumption that the bulge RGBB 
population will follow the same relations with metallicity as the Galactic GC system. We 
will show this to be a unique probe of the Galactic bulge stellar population. 
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Our comparison between the measured and predicted values for sightlines close to the 
plane are summarized in Table [TJ and that for the sightline (l,b) = (1,-6) in Table |2} We 
specify that while our calculation is done for a bulge MDF and compared to Galactic bulge 
observations, we expect the methodology to generalize well to future observations of the 
RGBB in kinematically-selected Gaia CMDs of the thin disk, thick disk and halo of the MW 



taken with Gaia observations (Lindegren 2010). 



We use the [Fe/H] and [a/Fe] of 204 bright RG stars toward Baade's window respectively 



measured by Zoccali et al. (2008) and Gonzalez et al. (2011a). The mean metallicity for these 



bright red giants is [M/H]= +0.047. However, for a composite stellar population such as the 
Galactic bulge, the metallicity distribution of one population (e.g. bright red giants) will 
not match the metallicity distribution of another population (e.g. RGBB stars), since the 
relative lifetimes are a function of metallicity. We correct for these two effects by using the 



relations from Renzini (1994): 



ALog t H B ' 
ALog t RG B 



-0.06ALogZ - 
-0.04ALogZ 



0.33AY, 



0.84AY, 



(20) 
(21) 



These relations predict that an HB star with [M/H]=0 and Y=0.27 will have a lifetime 
^hb ~17% longer than one with [M/H]= —1 and Y=0.25. For an RG star, the difference 
will be a reduction in lifetime of ~12%. In both cases the bulk of the difference comes from 
the metallicity component. For a standard helium to metals enrichment ratio of AY/AZ=1.5, 
the above equations reduce to: 



ALog tHB ~ +0.06ALogZ + 0.50AZ, 
ALog t RGB ~ -0.04ALogZ - 1.25AZ. 



(22) 
(23) 



The mean metallicity of RC stars has to be higher than that of the RG stars, since metal- 
rich RG stars have suppressed lifetimes whereas metal-rich HB stars have enhanced lifetimes, 
with the effect amplified by the monotonic relation between metallicity and initial helium 
abundance. A weighted mean must be computed: 



W it HB = exp 



(In 10)(0.06 + 0.04)[M/H] + (0.50 + 1.25)AZ 



LHB 



exp 



0.230[M/H] + 0.033(10 [M/H] - 1) 



(24) 



(25) 



The estimated mean metallicity of bulge RC stars is thus [M/H]= +0.074, which is expected 
to be a little higher than that of the bright red giants. 
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We then compute the mean metallicity distribution of the RGBB distribution by weight- 
ing over their respective relative lifetimes. For each RGBB star, we compute: 



W ijR GBB = /agbbI[m/h] * exp 



0.230[M/H] + 0.033(10 [M/H] - 1) 



(26) 



where the first factor of W^rgbb is the value of /J^b I [m/h] at that metallicity, the first-order 
approximation to the relative lifetime. The second factor corrects for the fact that metal-rich 
stars are numerically suppressed in the RG metallicity sample, as well as the fact that the 
normalization Nhb is enhanced at the high-metallicity end. The weighted-mean metallicity 
of the RGBB stars is thus [M/B] RGBB = +0.108. The predicted fraction is then: 



f RC I 

J RGBB I Bulge 



£ H B i 

/rgbbI[m/H]=o.io8 * exp 



0.0046 + 0.033(10 



(0.108-0.074) 



1) 



= (184/180). (27) 



The second factor corrects for the lower mean metallicity of the HB relative to that of the 
RGBB. The third factor, 184/180, accounts for the small number of BHB+RR+RHB stars 
that won't be included in a CMD selection box for the RC, estimated using HST proper 
motions toward the SWEEPS field (Clarkson et al. 2011). The term may be even larger 



if the bulge has an undiscovered extreme BHB population, a plausible outcome due to the 



UV-excess observed toward ellipticals and the bulges of disk galaxies (Terndrup & Walker 
1994). The derivation yields f B G BB \Bu\ge = (0.279 ± 0.015), if the assumption that the bulge 



and the Galactic GC system have the same input stellar physics holds. 

We also predict the brightness. As with the number counts, some empirical corrections 



are required due to the RC being a moving target. The RC is predicted by theory (Girardi & 



Salaris||2001 ) to decrease in brightness in J at a rate of 0.2 mag dex 1 , a prediction confirmed 



in observations of both the local Hipparcos population (Udalski 2000) and extragalactic 



systems (Pietrzyhski et al. 2010). 



We take the weighted mean of the brightness distribution of RGBB 

^ I RGBB I Bulge = W l>RGBB * {^gEbV + 0.2([M/H]|i - 0.074)), 



(2f 



where the first term gives the first-order weighted sum over the predicted brightness differ- 
ences, and the second term corrects for the second-order effect of the RC becoming fainter 
with increased metallicity, by bringing it back to its brightness at [M/H]= +0.074. We 
thereby obtain a mean brightness of AJ^ BB | B uige =0.767+0.038. 

However, we also estimate the impact of an additional factor. The RGBB-RC pair 



toward the bulge is fit using a double Gaussian (Nataf et al. 2011a Clarkson et al. 2011 



Nataf & Udalski 2011 Gonzalez et al. 2011b). Since the RGBB brightness distribution is 
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not Gaussian, some of the RGBB stars will be "ignored" by the fit due to the fact they are 
much closer in brightness to the RC than to the RGBB, this observational bias has the effect 
of lowering the observed number counts and increasing the observed brightness separation. 
We evaluated this effect in a multi-step simulation. We first use the observed metallicity 
distribution, and transform it into an RGBB brightness distribution, using the same weights 
as in the rest of this section. We add a Gaussian dispersion to each star to simulate the 
intrinsic width of the RGBB at a given metallicity using the relation of Section [81 as well as 



0.16 mag of Gaussian dispersion to account for the geometrical thickness of the bulge (Nataf 



& Udalski 2011) 



As can be seen in the bottom panel of Figure 17, the RGBB brightness distribution 
relative to the RC is predicted to be skewed. An actual measurement using the double- 
Gaussian method will be biased, by systematically missing out on the RGBB stars that 
are brighter than or near the brightness of the RC. We test for this observational bias in 
two ways, the first by simply adding the RC, and the second by adding the RC and a RG 



luminosity function. In the first case (third panel of Figure 17), AI B q BB is shifted to 0.837 
mag, and Jrgbb drops to 0.241. The predicted dispersion for the RGBB is 0.258 mag. In 



the second case (bottom panel of Figure 17), AI§q BB is shifted to 0.844 mag, and f^GBB 
shifts to 0.259. The predicted dispersion for the RGBB is 0.280 mag. Adding an RC shifts 
the measured RGBB peak to a fainter luminosity and reduces its number counts. Adding an 
RG luminosity function then mostly restores the RGBB number counts. This is due to the 
fact that a symmetric Gaussian is being fit to a skewed RGBB. The Gaussian then "scoops" 
up some of the RG stars as a byproduct of enforcing it symmetric profile. 



These parameters are inconsistent with those measured in Nataf et al. (2011a) and 



revised in Section 10 The brightness peak is brighter by 0.1 mag, the brightness dispersion 
is lower, and the number counts are smaller by 20%. This presents a strong case that the 
input physics for the bulge stellar population (age, helium, etc) are different from those of 
the Galactic GC system. The measured and simulated luminosity function for the Galactic 
bulge, with their stark differences, are shown in Figure [17} 
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Table 1: Observable parameters for the RGBB toward a triaxial ellipsoid sightline of 
the Galactic bulge. The second column lists the measured values. The third column 
predicts the values of the brightness, normalization and brightness dispersion for the RGBB. 
The fourth column predicts what values would be measured by the double-Gaussian method. 



Parameter 


Measured 


Predicted 


Predicted with RC+RG+RGBB LF 


A T RC 


0.737±0.012 


0.767±0.038 


0.844±0.043 


fHB 
J RGBB 


0.201±0.012 


0.279±0.015 


0.259±0.017 


&RGBB 


0.220±0.010 


0.279±0.015 


0.280±0.015 



Table 2: Observable parameters for the RGBB toward the double-RC sightline 
(l,b) = (1,-6). The second column predicts the values of the brightness, normaliza- 
tion and brightness dispersion for the RGBB. The third column predicts what values would 
be measured by the double-Gaussian method. 



Parameter 


Predicted 


Predicted with RC+RG+RGBB LF 


A jRC 
RGBB 


0.658±0.034 


0.745+0.038 


fHB 
J RGBB 


0.265±0.016 


0.256+0.015 


CRGBB 


0.281±0.015 


0.234+0.012 
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-1 

^gbb-^c Predicted 



Fig. 17. — From top. Panel 1: The metallicity distribution toward Baade's window from 



Gonzalez et al. (2011a). Panel 2: Predicted brightness distribution for the RGBB relative 



to the mean of the RC in /. Panel 3: Combined predicted brightness distribution for the 
RGBB and RC. Panel 4: Predicted brightness distribution with a RG luminosity function. 
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Fig. 18. — TOP: Measured luminosity function for the Galactic bulge RC+RGBB+RG 
branch. BOTTOM: Predicted luminosity function for the Galactic bulge RC+RGBB+RG 
branch given the bulge metallicity distribution and the metallicity relations measured in 
Galactic GCs. 



We repeat the exercise for the metallicity distribution toward (l,b) = (1, —6). (Zoccali 



et al. 2008 Gonzalez et al. 2011a), which has a mean metallicity in bright red giants of 



[M/H]= —0.067. This sightline is important due to the recent discovery that the Galactic 



bulge RC bifurcates at large latitudes, likely due to an X-shaped bulge (Nataf et al. 2010 
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McWilliam & Zoccali 2010). The RGBB toward these sightlines must be properly understood 



if there is to be any prospect of accurately fitting for the degenerate RCs. 
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Fig. 19. — Same as Figure 17 but corresponding to the metallicity measurements toward the 
sightline (l,b) = (1,-6). 



We simulate these sightlines as having negligible geometrical dispersion and an intrinsic 



brightness dispersion for the RC of 0.17 mag in / (Nataf & Udalski 2011). The mean 



brightness of the RGBB stars is 0.658 mag fainter than the RC, with fjiQ BB = 0.265. After an 
RG+RC luminosity function is added, AI B q BB = 0.745, ctrgbb = 0.234 mag, and f B Q BB = 
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0.256. If the Galactic bulge RGBB indeed has over 20% of the numbers counts of the 
Galactic bulge RC within the Milky Way's X-wings, then it must be taken into account 
when modelling those sightlines. In particular, it could explain why there is a low measured 
difference in radial velocity distributions between the two RCs (De Propris et al. 2011): 
perhaps the brighter RGBB is heavily mixed with the fainter RC. 



11.1. RGBB Properties for the Bulge Solution: Enhanced Helium 
Enrichment, Possibly with Some Age Variation 

We run a few illustrative stellar models using the Yale Rotating Evolution Code (YREC) 



with diffusion (Sills et al. 2000 Delahaye et al. 2010|) and empirically calibrated bolometric 

a = 1.922 



corrections ( Pinsonneault et al. 2004 An et al.]|2007 ). The mixing length is set to 
to match current data of the solar radius, luminosity, and atmospheric metals to hydrogen 



ratio (Grevesse & Sauval 1998). We test two hypotheses that might explain the increased 



brightness and decreased number counts of the Galactic bulge RGBB relative to the Galactic 
GC calibration. The first is that the Galactic bulge has enhanced helium enrichment without 
a younger age, and the second is that the Galactic bulge has a younger age with standard 
helium enrichment. The results are summarized in Table [3] and shown in Figure 20 



Table 3: Observational discrepancy compared to predicted evolutionary effects of enhanced 
enrichment helium and a younger age. 



Parameter 


Observational Discrepancy 


AY = +0.06 


At = -5 Gyrs 


A tRC 
RGBB 


-0.107±0.045 


-0.11 


-0.29 


fHB 
J RGBB 


-22.4±8.0% 


-32.2% 


-25.6% 



Both adjustments do an effective job of matching the observational discrepancy on 
number counts. However, the helium-enhanced model yields a superior fit to the brightness 
measurement: adding AY= +0.06 makes the RGBB brighter by 0.11 mag, whereas reducing 
the age by ~5 Gyr increases the brightness by 0.29 mag. The observational discrepancy of 
0.107+0.045 mag is much more well-matched by the former. 

Ultimately, both factors may play a role. A complete analysis of the anomalous RGBB 
would require a suite of stellar tracks across the composite metallicity range of the Galactic 
bulge, as well as an accounting of the fact that the RC lifetime and brightness will also vary 
with the age and initial helium abundance of a stellar population - though not as steeply as 
for the RGBB. 
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Fig. 20. — We compare a canonical bulge stellar track with [M/H]= +0.20, Y=0.29 and 
t=12 Gyrs, shown in red with a helium-enhanced track of the same age and metallicity in 
the top left panel, and a younger track of the same helium and metallicity in the bottom 
left panel, the latter two shown in blue. The right panels show the corresponding luminosity 
functions. Both modifications yield a ~30% reduction in the lifetime of the RGBB, but the 
age modification yields a change in brightness 3x that obtained by increasing the helium. 
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Since there are three measurable parameters for the Galactic bulge RGBB, the prospects 
of tightly constraining both the age-metallicity and helium-metallicity relations, to yield a 
Galactic bulge age-helium-metallicity relation, are decent. The prospects improve further if 
the RGBB properties are found to be measurable toward multiple sightlines of the Galactic 
bulge that differ in their metallicity distributions. 



11.2. Predicted RGBB Properties for the Bulge: Caveats 



However scientifically satisfactory it may be to have the GC relations with metallicity 
for the RGBB-RHB pair with which to to construct a null hypothesis to test for the bulge, 
it must be pointed out that this anchor is itself imperfect, as it is not precisely known. A 
significant scientific concern is the underlying foundation of this approach - the assumption 
that there are pertinent mean relations for the GCs. These are diverse stellar systems with 



significant variations in self-enrichment profiles ( 


Carrettaet al.|2010 


Caloi&D'Antona|2011 


Valcarce & Catelan 


2011 


), in age ( 


Marm-Franch et al. 


2009 


Dotter et al. 


2011 


), and other 



properties. Milone et al. (2011) used observations of 47 Tuc in 9 bandpasses to estimate 
that ~70% of the stars are helium and nitrogen enriched as well as oxygen deficient. If this 
is found to be the norm for the Galactic GC system, it will be necessary to incorporate 
the expectation that the RGBB in Galactic GCs should be a little brighter, and a little 
less-populated, than from a population with canonical abundance patterns. 

The weights used are not a significant source of systematic error. We repeat the same 



calculation with the nearly pure RC spectroscopic sample measured by Hill et al. (2011), 
which does not require weights derived from stellar evolution models as it is already anchored 
at the core helium-burning phase. We find that AI B q BB increases by 0.040 mag, and J B gbb 
increases by 0.006. The source of this difference is that the sample of Hill et al. (2011) 



report slightly higher metallicities, which they speculate in their Section 4.1.4 may be due to 
upgrades in their spectroscopic reduction procedure. In this paper we have used the combined 



sample of Zoccali et al. (2008) and Gonzalez et al. (2011a) to have a uniform sample for the 



sightlines toward the triaxial ellipsoid component of the bulge and the X-shaped component 
of the bulge. 



The predicted value of /J&§b_b| Bulge may be artificially decreased by systematic effects. 



Hill et al. (2011) argue that spectroscopy of bulge giants may systematically underestimate 



the metallicities of the most metal-rich stars. If true, this would increase the predicted values 
of AIftQ BB and f^GBB j increasing the discrepancy with the observed values. 



We also point out that the metallicity distribution of the bulge is turning out to be 
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more complex than previously assumed, as it is bimodal, and possibly more complex. |Hill 



et al. (2011) found two peaks in their sample, one at [Fe/H]= —0.30 and the other at 
Fe/H]= +0.32. When we use their MDF, we do not find a significant change in urgbb- 



Bensby et al. (2011), using observations of microlensed Galactic bulge dwarf and subgiant 



stars, also find two peaks, but that are more broadly separated than those measured from 
RC stars, at [Fe/H]= —0.60 and [Fe/H]= +0.30. They also show that fits of their measured 
spectroscopic temperatures and gravities to isochrones imply a significant spread in age. We 
expect that the calculations performed in this section will likely need to be repeated in the 
future as observational constraints on the bulge's age and metallicity continuously improve 
with time. 



12. Discussion and Conclusion 

In this work, we have introduced and justified a scientifically robust parametrization 
with which to study the RGBB and the associated RG luminosity function. The relevant pa- 
rameters can be fit concurrently rather than sequentially, and the use of maximum-likelihood 
estimation removes any concern that the size and position of bins could distort the output 
parameters. By combining this parametrization with the photometry from the ACS GC 
survey ( |Sarajedini et~aL][2007| and that of the WFPC2 GC survey QPiotto et aL||2002| ), we fit 



for the brightness and color of the RGBB in 72 GCs, the brightness and color of the MSTO 
in 55 GCs, and the brightness of the RHB in 31 GCs. We also fit for the strength of the 
RGBB, EWrgbb and the number of HB stars for all 72 GCs. There are several empirical 
achievements in this work. 

This is the most robust investigation of RGBB star counts in GCs. Measurements of 
Nrgbb reach precisions of 10% in the most populous clusters, and the mean relations for 
EWrgbb and /rgbb are determined to 4% accuracy. Previous investigations had relied on 
the Reump parameter, a composite measure of the RG luminosity function and the strength 
of the RGBB, that has a lower signal to noise ratio. The RGBB in 47 Tuc, previously 
detected with ~5cr-significance, is now detected with ~9a-significance due to this different 
parametrization. Measurements of the strength of the RGBB feature are now on firm-enough 
footing that it is in itself a tool with which to precisely compare GC observations to stellar 
model predictions. 

We also compute predicted values of AI§q BB and /rgbb f° r two Galactic bulge sight- 
lines given the assumption that these parameters have the same functional dependence on 
metallicity in the bulge as they do for GCs. The results are not consistent with those found 



in Nataf et al. (2011a) and revised in this work - the predicted RGBB luminosity is fainter, 
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broader in magnitude spread, and more significant in number counts. As discussed in |Nataf 



et al. (2011a), one path to resolve this discrepancy is to posit enhanced helium-enrichment 



for the Galactic bulge. A higher value of AY/AZ would make the RGBB stars brighter, 
thereby decreasing the size of the derivative of brightness with metallicity, and it would also 
decrease the lifetime of the RGBB. Enhanced helium for the bulge is also a prediction of 



chemical evolution, due to the enhanced a-element abundance toward the bulge (Catelan 



2007). The question of exactly how much helium is needed will be tackled in a future paper, 



where we will compare these results to stellar model predictions. 

Our analysis of number counts may lead to a resolution of a longstanding issue in RGBB 
astrophysics, that the observed brightness of the RGBB in GCs is ~0.2 mag fainter than 



that predicted by models (Fusi Pecci et al. 1990 Di Cecco et al. 2010 Cassisi et al. 2011 



Troisi et al. 2011). Including overshooting beyond the formal boundary of the convective 



envelope has been proposed as a solution (Alongi et al. 1991). If this is modification to 



stellar models yields a better match to data due to genuine processes in stars and not due to 
a coincidence, then there is the straightforward prediction that stellar models with adjusted 
overshoot should also yield a better match to number counts as well as the brightness peak 



of the RGBB. Separately, Bjork & Chaboyer (2006) argued that the theoretical uncertainties 



in the mixing length, low-temperature opacities and the observational uncertainties in the 
abundance of a-capture elements could adjust the predicted number counts of the RGBB to 
no more than crRBump ~0.01, equivalent to a ~3% uncertainty in the lifetime of the RGBB. 
These predictions of stellar evolution are now testable. 

In the next decade, large-scale surveys and improved-instrumentation will further the 
depth of astrophysical research accessible with the RGBB. Observations of this galaxy by 
the Gaia mission^ will allow detailed investigations of the RGBB in the disk of this galaxy. 
Meanwhile, if 30 meter telescopes are built, higher-quality CMDs of GCs will be available 
throughout the local group. Though we do not expect split main-sequences to be as observ- 
able as uj Cen is using HST (Bedin et al. 2004) , a split RGBB should be observable toward 



those kinds of systems in Andromeda or Triangulum if they exist. Indeed, the split RGBB 



of u Cen can be viewed from ground-based, 1-meter telescopes without adaptive optics (Rey 



et al.|[2004 ). As much as the astrophysics accessible with HST observations of the RGBB in 
GCs is an upgrade over what was previously available, we forsee even more significant gains 
in the coming decade due to observational efforts listed above. 
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A. Testing our Methodology with Monte Carlo Methods 

We use Monte Carlo simulations to determine the statistical robustness of our method- 
ology. We compare the relative diagnostic power of the parameter EWrgbb and RBump, and 
we investigate whether or not out maximum-likelihood approach gives reliable estimates of 
the luminosity function parameters and their errors. 

We sample from a population of stars drawn from the probability density function 
characterized by Equation We set B = 0.72, EWrgbb = 0.30, cr RGBB = 0.05, and 
ttirgbb = 0, where we use the generalized notation kirgbb for this section, and not Irgbb 
or Vrgbb as used elsewhere in this paper, as the RGBB could in principle be investigated 
with other bandpasses. The stars are distributed within the range ttirgbb — 1-5 < m < 
wirgbb + 7.0, with the total number of stars distributed log-uniformly between ~3000 and 
~300000. We keep the ~3% of the stars within the range tjirgbb — 1.5 < m < ttirgbb + 2.0. 
We then run an MCMC on the simulated luminosity function in the same manner as in 
Section [3j We do so with the same observatiional biases, keeping only the outputs that have 
both a measured an expected value of N RG bb > 3.0. Over the 3,000 runs of the simulation, 
half have their MCMCs run without priors, and the other half have the same priors as those 
we imposed on the silver sample. 

The larger initial range is used so that bins in the selected range have uncorrelated 
rather than anti-correlated number counts. In a real stellar population, the number of stars 
in a magnitude interval (mi, mi + 5m) is independent of the number of stars in another 
interval (m2, wi2 + 8m) - both populations are distributed as their own independent Poisson 
random variables. This will not be true in a simulation with a fixed total number of stars, 
where 1 additional star in 1 bin means 1 less available star for the other bins. We are 
effectively breaking this anti-correlation by sampling from a population reservoir that is ~30 
times more populated. 
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A.l. The Relative Diagnostic Power of the EWrgbb an d RBump Parameters 

We find that the parameter EWrgbb averages approximately twice the statistical sig- 
nificance of the Reurnp parameter. The behavior of the two parameters is shown in Figure 

That the ratio of significances would be near 2 is not surprising. We demonstrated 



in Section 32 that for a typical RG+RGBB luminosity function, ~73% of the stars in the 
interval Vrgbb — 0.5 < V < Vrgbb + are RG stars and not RGBB stars. Since statistical 
fluctuations in the number of stars scales as the square root of the number of stars, the 
factor of ~2 decrease in statistical significance for RBump is in fact the expected outcome of 
effectively diluting the sample size by a factor of 4. Other sources of error, such as that of 
the denominator in the RBump parameter, are also present but they are not dominant. 

The behavior of RBump at low number counts is particularly devastating. The statistical 
significance is frequently below zero, which would imply an unphysical negative normalization 
for the RGBB. The significance only reaches 1 in the 50th percentile when Nrgbb surpasses 
10. 
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Fig. 21. — The statistical significance of the RGBB given the choice of parametrization, as 
a function of the true number of RGBB stars in the model distribution. Filled black squares 
denote the median statistical significance for EWrgbb with the error bars denote 32nd and 
68th percentiles. Similarly for the empty black squares and Reump- 



A. 2. The Inferred Population Parameters Versus the True Population 

Parameters 

We find that the maximum likelihood approach yields an unbiased estimator for EWrgbb 
at large number counts, but one which is biased toward higher inferred values of EWrgbb 
for Nrgbb ^ 10. The bias disappears if we impose the same priors that we imposed on 
the silver sample, which was constructed out of GCs that had a measured best-fit value 
N RG bb < 10. 

The maximum-likelohood value of uirgbb is found to be an unbiased estimator of the 
true value of wirgbb- However, the scatter increases for lower values of Nrgbb- At low 
number counts, any methodology is at risk of fitting to other peaks in the distribution that 
arise from statistical fluctuation. The lack of bias is due to the fact that these other peaks 
need not be either fainter or brighter than the true peak of the RGBB. 
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Fig. 22. — TOP: The distribution of differences between the value of EWrgbb inferred by 
the MCMC, and the true value from which the histogram is constructed. Error bars denote 
32nd and 68th percentile. Empty circles denote the distribution of inferences without priors, 
filled circles with priors. BOTTOM: Same as top, but for the parameter iurgbb- 

Both parameter comparisons demonstrate the urgent need for a broad investigation of 
the RG+RGBB luminosity function over the full pertinent range of ages, metallicities and 
helium enrichments, to ascertain the theoretical robustness of these priors. 



A. 3. The Inferred Errors in the Population Parameters Versus the True 
Errors in the Inferred Population Parameters 

It is important to demonstrate not just that our parameter estimates are unbiased 
in the mean, but that the errors in our parameter estimates are unbiased as well. The 
standard deviation of the differencce between the maximum-likelihood value and the true 
value of the parameters should be equal to mean of the errors reported. We find that 
reliable determinations of the errors are obtained by our maximum-likelihood method for 
Nrgbb > 10, regardless of the use of priors. At low number counts, the errors in EWrqbb 
remain unbiased with the priors we used to construct our silver sample, but a small bias 



remains in the errors in the inferred brightness. The results are shown in Figure 23 
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Fig. 23. — TOP LEFT: The 32nd, 50th, and 68th percentile of the error in EWrgbb inferred, 
for MCMCs with priors (filled circles) and without (empty circles). BOTTOM LEFT: The 
standard deviation of the difference between the true EWrgbb and the value measured by 
the maximum-likelihood method. TOP RIGHT: The 32nd, 50th, and 68th percentile of the 
error in ttirgbb inferred, for MCMCs with priors (filled circles) and without (empty circles). 
BOTTOM RIGHT: The standard deviation of the difference between the true m RGBB and 
the value measured by the maximum-likelihood method. 

There is a simple explanation for this behavior. At low number counts, the MCMC risks 
jumping from the true RGBB brightness peak to other statistical fluctuations that may crop 
up. The measured errors will then be the errors in the position and normalization of that 
peak, rather than of the true peak. 

Whereas the best-fit values of EWrgbb and turgbb are unbiased at low number counts, 
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the error in muGBB is likely underestimated by the MCMC for N RGBB < 10. The error 
reported is the error in the peak which is fit for, and not the difference between the location 
of the true peak and that of the peak which is fit for. 

In practice, the relation for the brightness of the RGBB with metallicity will allow 
astronomers to rule out peaks that differ from the true peak by 1.0 mag or more. However, 
this is not possible in the less frequent case where a secondary peak shows up within 0.1 
mag or less. Additionally, other catastrophic failures of fitting at low number counts, such 
as peaks with grgbb > 0.25, or EWrgbb > 1-0, will be selected against. 
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Table 4: Cluster metallicities, estimated V-band apparent distance modulus, Vrgbb, as 
well as (V - I) rgbb and (B — V)rgbb where measured. Clusters with an asterix in their 
name are part of the silver sample, l^-band magnitudes have for GCs from the WFPC2 
dataset have been made brighter by 0.0365 mag, in the manner described in Section [2] The 
(B — V) measurements are not adjusted. 



Name 


[Fe/H 




[M/HJ 


(m-Mj V 


t r 

'RGBB 


(V - I) 


S D T/"\ 

[B - V) 


ARP0002* 


-1.74±0.08 


-1.47±0.10 


17.59 


17.940±0.044 


1.08 




IC04499* 


-1.62±0.09 


-1.36±0.11 


17.08 


17.332±0.018 


1.24 




LYNGA07 


-0.62±0.10 


-0.38±0.11 


16.78 


18.039±0.017 


1.93 




NGC 104 


-0.76±0.02 


-0.45±0.05 


13.37 


14.507±0.005 


1.03 




NGC 1261 


-1.27±0.08 


-1.02±0.10 


16.09 


16.599±0.010 


0.97 




NGC 1851 


-1.18±0.08 


-0.90±0.10 


15.47 


16.078±0.009 


1.00 




NGC 1904* 


-1.58±0.02 


-1.36±0.06 


15.59 


15.877±0.017 




0.82 


NGC 2808 


-1.18±0.04 


-0.94±0.07 


15.59 


16.235±0.013 


1.21 




NGC 3201 


-1.51±0.02 


-1.27±0.05 


14.20 


14.649±0.032 


1.30 




NGC 362 


-1.30±0.04 


-1.09±0.07 


14.83 


15.399±0.007 


0.98 




NGC 4590* 


-2.27±0.04 


-2.02±0.09 


15.21 


15.149±0.011 


1.06 




NGC 4833* 


-1.89±0.05 


-1.62±0.07 


15.08 


15.246±0.029 


1.39 




NGC 5024 


-2.06±0.09 


-1.79±0.11 


16.32 


16.488±0.020 


0.97 




NGC 5272 


-1.50±0.05 


-1.26±0.07 


15.07 


15.448±0.010 


0.96 




NGC 5286 


-1.70±0.07 


-1.43±0.09 


16.08 


16.287±0.009 


1.28 




NGC 5634 


-1.93±0.09 


-1.66±0.11 


17.16 


17.371±0.030 




0.78 


NGC 5824 


-1.94±0.14 


-1.67±0.15 


17.94 


18.084±0.028 




0.88 


NGC 5904 


-1.33±0.02 


-1.05±0.05 


14.46 


14.963±0.009 


1.01 




NGC 5927 


-0.29±0.07 


-0.06±0.09 


15.82 


17.233±0.014 


1.55 




NGC 5986 


-1.63±0.08 


-1.37±0.10 


15.96 


16.397±0.030 


1.32 




NGC 6093 


-1.75±0.08 


-1.58±0.10 


15.56 


15.999±0.020 


1.24 




NGC 6101* 


-1.98±0.07 


-1.71±0.09 


16.10 


16.235±0.021 


1.11 




NGC 6139 


-1.71±0.09 


-1.44±0.11 


17.35 


17.867±0.019 




1.49 


NGC 6144* 


-1.82±0.05 


-1.55±0.07 


15.86 


16.099±0.030 


1.51 




NGC 6171 


-1.03±0.02 


-0.66±0.07 


15.05 


15.870±0.038 


1.50 




NGC 6205 


-1.58±0.04 


-1.36±0.07 


14.33 


14.774±0.013 


0.98 




NGC 6218 


-1.33±0.02 


-1.03±0.06 


14.01 


14.785±0.011 


1.21 




NGC 6229 


-1.43±0.09 


-1.17±0.11 


17.45 


17.899±0.025 




0.84 


NGC 6235* 


-1.38±0.07 


-1.12±0.09 


16.26 


16.763±0.014 




1.10 
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NGC 6254 


-1.57±0.02 


-1.30±0.05 


14.08 


14.787±0.012 


1.31 




NGC 6284 


-1.31±0.09 


-1.06±0.10 


16.79 


17.370±0.050 




1.10 


NGC 6304 


-0.37±0.07 


-0.14±0.09 


15.52 


16.904±0.014 


1.62 




NGC 6316* 


-0.36±0.14 


-0.13±0.15 


16.77 


18.181±0.046 




1.47 


NGC 6341 


-2.35±0.05 


-2.01±0.07 


14.65 


14.666±0.013 


0.99 




NGC 6352 


-0.62±0.05 


-0.48±0.07 


14.43 


15.732±0.022 


1.33 




NGC 6356 


-0.35±0.14 


-0.12±0.15 


16.76 


18.076±0.016 




1.17 


NGC 6362 


-1.07±0.05 


-0.82±0.07 


14.68 


15.485±0.021 


1.06 




NGC 6366* 


-0.59±0.08 


-0.35±0.10 


14.94 


16.077±0.019 


1.90 




NGC 6388 


-0.45±0.04 


-0.30±0.07 


16.13 


17.650±0.010 


1.45 




NGC 6397* 


-1.99±0.02 


-1.73±0.06 


12.37 


12.533±0.046 


1.18 




NGC 6402 


-1.39±0.09 


-1.13±0.11 


16.69 


17.317±0.025 




1.42 


NGC 6426* 


-2.26±0.10 


-1.98±0.11 


17.68 


17.736±0.019 


1.44 




NGC 6440* 


-0.20±0.14 


0.03±0.15 


17.95 


19.431±0.021 




1.97 


NGC 6441 


-0.44±0.07 


-0.29±0.09 


16.78 


18.395±0.008 


1.60 




NGC 6496* 


-0.46±0.07 


-0.23±0.09 


15.74 


16.975±0.015 


1.30 




NGC 6539* 


-0.53±0.14 


-0.21±0.15 


17.62 


18.847±0.024 




1.81 


NGC 6541 


-1.82±0.08 


-1.50±0.10 


14.82 


15.029±0.014 


1.10 




NGC 6569 


-0.72±0.14 


-0.48±0.15 


16.83 


17.781±0.019 




1.39 


NGC 6584 


-1.50±0.09 


-1.24±0.11 


15.96 


16.342±0.013 


1.07 




NGC 6624 


-0.42±0.07 


-0.19±0.09 


15.36 


16.617±0.013 


1.34 




NGC 6637 


-0.59±0.07 


-0.37±0.09 


15.28 


16.394±0.013 


1.20 




NGC 6638 


-0.99±0.07 


-0.74±0.09 


16.14 


17.038±0.015 




1.23 


NGC 6642* 


-1.19±0.14 


-0.94±0.15 


15.79 


16.560±0.023 




1.23 


NGC 6652 


-0.76±0.14 


-0.52±0.15 


15.28 


16.366±0.032 


1.15 




NGC 6656* 


-1.70±0.08 


-1.42±0.10 


13.60 


13.974±0.013 


1.42 




NGC 6681 


-1.62±0.08 


-1.36±0.10 


14.99 


15.627±0.019 


1.08 




NGC 6717* 


-1.26±0.07 


-1.01±0.09 


14.94 


15.729±0.021 


1.23 




NGC 6723 


-1.10±0.07 


-0.72±0.09 


14.84 


15.611±0.009 


1.07 




NGC 6752 


-1.55±0.01 


-1.23±0.05 


13.13 


13.625±0.020 


1.03 




NGC 6760 


-0.40±0.14 


-0.17±0.15 


16.72 


18.284±0.020 




1.66 


NGC 6809* 


-1.93±0.02 


-1.62±0.06 


13.89 


14.151±0.013 


1.10 




NGC 6838* 


-0.82±0.02 


-0.53±0.05 


13.80 


14.833±0.018 


1.30 




NGC 6864 


-1.29±0.14 


-1.04±0.15 


17.09 


17.678±0.011 




1.00 



NGC 6934 -1.56±0.09 -1.30±0.11 16.28 16.648±0.013 1.09 



NGC 6981 



-1.48±0.07 -1.22±0.09 16.31 



16.715±0.017 1.02 
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NGC 7006 


-1.46±0.06 


-1.20±0.08 


18.23 


18.641±0.004 


1.06 




NGC 7078 


-2.33±0.02 


-2.04±0.08 


15.39 


15.315±0.021 


1.06 




NGC 7089 


-1.66±0.07 


-1.36±0.09 


15.50 


15.815±0.008 


1.02 




NGC 7099* 


-2.33±0.02 


-2.06±0.07 


14.64 


14.712±0.020 


1.01 




PYXIS00* 


-1.40±0.10 


-1.14±0.11 


18.63 


19.156±0.032 


1.35 




RUPR106* 


-1.78±0.08 


-1.51±0.10 


17.25 


17.489±0.042 


1.19 




TERZAN8* 


-2.00±0.20 


-1.73±0.21 


17.47 


17.660±0.019 


1.12 
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Table 5: Apparent Brightness and Color for the Main-Sequence Turnoffs for all clusters 
observed in the ACS survey, and relative brightness of the RHB for all clusters with a 
prominent RHB. Clusters with an asterix in their name are part of the silver sample. V- 
band magnitudes have for GCs from the WFPC2 dataset have been shifted in the manner 
described in Section [2] 



Name 


' M o 1 U 


V y ± ) IVlblU 


4-1 y RGBB 


A jRHB 

RGBB 


A t? Pnnno* 
Arvr UUUZ 


01 KQQ 
Zl.OOO 


n 7nn 
U. lUO 


o.oyo±U.U44 




TPn/l A OO* 


01 ns7 
Zl.Uo ( 


n fic;o 
U.ooy 


o. I 00±U.Ulo 






on ofio 
zu.ysz 


i fi/i i 
1.041 


o oa Q-un m 7 
z.y4o±u.ui / 


n c^no-un m o 
U.oUZ±U.Uiy 


MP p in/1 
INLrU 1U4 


1 < . < 1U 


n soft 
u.oyo 


q onQ-un nnc; 
O.ZUoitU.UUO 


n qfi/i-un nnfi 
U.oo4±U.UU0 


MP p i ofil 
INLjU IZOl 


on lie; 
ZU. 110 


n kq /i 


QKi fi_i_n m n 
o.olo±U.UlU 


n q 1 a 4~n m 1 
-U.ol4±U.Ul 1 


MP p 1QC1 

INLjU loOl 


1 n c;fiq 
ly.ooo 


U.OOO 


q a ocin nnn 
o.4oo±U.UUy 




mpp osns 

INLrU ZOUo 


i n fioi 
ly.oy 1 


n Qc;n 
U.ooU 


q a c;fi_i_n ni q 
o.400±U.Ulo 


n ons-i~n m q 
-U.zUoztU.Ulo 


mpp Qom 
INLjU ozUl 


1 S OOfi 

lo.zzo 


n n/i i 
U.y41 


q K77xn nqo 
o.O 1 1 ztU.Uoz 


n on/i -i~n c\a 7 
-U.zy4±U.U4 1 


MPP qfiO 
IN lor U OOZ 


1 Q OA 1 

lo.o41 


n fins 
U.oUo 


q a a q~i~n nn7 
o.44o±U.UU 1 


n oc;fi-i~n nna 
-U.zooztU.UUo 


mp p a c;on* 


iy.Uo4 


U.oyo 


o.yoOztU.Ull 




MPP /lOQQ* 

INlorO 4ooo 


iy. loo 


1 noo 
l.UZZ 


q cno-un non 
o.oyz±u.uzy 




MPP cnO/1 
iNlorlo 0UZ4 


on Q1 Q 
ZU.olo 


n c;7o 
U.O 1 z 


q oo^-un non 
o.oZ0±U.UZU 




MPP K070 

InIotLv ozv z 


i n n7o 
ly.U ( y 


n K7n 
U.o < y 


q fiqi _i_n m n 
o.OolihU.UlU 




MPP KOOfi 
InIotIo DZOO 


on 1 fin 
zU. loU 


n nQ 
U.yoO 


q 2704-0 nnn 
o.o i z±U.UUy 




MPP cnn/1 
InoU 0yU4 


lo.4 i 4 


n fiOQ 
U.ozo 


q ki o-un nnn 
o.olz±U.UUy 


n qfinj~n ni fi 
-U.ooUztU.Ulo 


MPP Cn07 

invjU oyz i 


on 1 Q7 
ZU. lo ( 


1 ooo 
i.zzy 


o on/i -i~n m a 
Z.yU4±U.U14 


U.Oo 1 ±U.U10 


NGC 5986 


20.160 


0.980 


3.763±0.030 




NGC 6093 


19.863 


0.876 


3.864±0.020 




NGC 6101* 


20.089 


0.709 


3.854±0.021 




NGC 6144* 


19.957 


1.176 


3.858±0.030 




NGC 6171 


19.342 


1.190 


3.472±0.038 


0.025±0.040 


NGC 6205 


18.529 


0.597 


3.755±0.013 




NGC 6218 


18.328 


0.871 


3.543±0.011 




NGC 6229 








-0.320±0.031 


NGC 6254 


18.538 


0.932 


3.751±0.012 




NGC 6304 


19.878 


1.333 


2.975±0.014 


0.586±0.014 


NGC 6316* 








0.285±0.046 


NGC 6341 


18.575 


0.567 


3.909±0.013 




NGC 6352 


18.805 


1.011 


3.073±0.022 


0.450±0.023 
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NGC 6356 








0.531±0.017 


NGC 6362 


18.896 


0.719 


3.412±0.021 


-0.052±0.024 


NGC 6366* 


19.081 


1.602 


3.004±0.019 


0.435±0.029 


NGC 6388 


20.797 


1.143 


3.147±0.010 


0.309±0.013 


NGC 6397* 


16.546 


0.799 


4.013±0.046 




NGC 6426* 


21.695 


1.073 


3.959±0.019 




NGC 6440* 








0.622±0.022 


NGC 6441 


21.511 


1.318 


3.116±0.008 


0.418±0.011 


NGC 6496* 


19.933 


0.985 


2.957±0.015 


0.471±0.016 


NGC 6539* 








0.418±0.026 


NGC 6541 


18.814 


0.735 


3.785±0.014 




NGC 6569 








0.141±0.021 


NGC 6584 


20.027 


0.686 


3.686±0.013 




NGC 6624 


19.713 


1.018 


3.096±0.013 


0.460±0.014 


NGC 6637 


19.561 


0.881 


3.167±0.013 


0.338±0.013 


NGC 6638 








0.046±0.019 


NGC 6652 


19.512 


0.818 


3.146±0.032 


0.321±0.033 


NGC 6656* 


17.854 


1.069 


3.880±0.013 




NGC 6681 


19.230 


0.716 


3.603±0.019 




NGC 6717* 


19.308 


0.901 


3.579±0.021 




NGC 6723 


19.072 


0.724 


3.461±0.009 


-0.078±0.014 


NGC 6752 


17.384 


0.654 


3.758±0.020 




NGC 6760 








0.479±0.021 


NGC 6809* 


17.939 


0.716 


3.789±0.013 




NGC 6838* 


17.971 


0.964 


3.138±0.018 


0.295±0.019 


NGC 6864 








-0.117±0.013 


NGC 6934 


20.276 


0.711 


3.629±0.013 




NGC 6981 


20.320 


0.640 


3.604±0.017 




NGC 7006 


22.262 


0.685 


3.621±0.004 




NGC 7078 


19.269 


0.651 


3.954±0.021 




NGC 7089 


19.500 


0.639 


3.685±0.008 




NGC 7099* 


18.658 


0.597 


3.946±0.020 




PYXIS00* 


22.692 


0.956 


3.537±0.032 


-0.309±0.036 


RUPR106* 


21.105 


0.798 


3.616±0.042 


-0.615±0.048 


TERZAN8* 


21.554 


0.721 


3.894±0.019 
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Table 6: Number density parameters for the RGBB are reported, as well as raw number 
counts for the RGBB and the HB. Clusters with an asterix in their name are part of the 
silver sample. 



AT 

JName 


J~l TI T 

&VV RGBB 


AT 

1N| RGBB 


AT 

Nhb 


fHB 
J RGBB 


ARP0002* 


0.370 ± 0.209 


6.0 ±3.3 


25 


0.239 ±0.139 


IC04499* 


0.110 ±0.071 


4.6 ±2.9 


103 


0.045 ± 0.029 


LYNGA07 


0.412 ±0.104 


30.3 ±6.8 


59 


0.517 ±0.129 


NGC 104 


0.322 ±0.041 


122.3 ± 14.2 


545 


0.224 ± 0.028 


NGC 1261 


0.213 ±0.063 


22.5 ±6.3 


231 


0.097 ±0.028 


NGC 1851 


0.294 ±0.063 


51.7 ±9.7 


397 


0.130 ±0.025 


NGC 1904* 


0.155 ±0.069 


10.0 ±4.3 


168 


0.059 ± 0.026 


NGC 2808 


0.347 ±0.054 


170.3 ±22.8 


1325 


0.129 ±0.018 


NGC 3201 


0.477 ±0.183 


14.9 ±5.1 


69 


0.215 ±0.078 


NGC 362 


0.306 ±0.065 


44.8 ± 8.4 


368 


0.122 ±0.024 


NGC 4590* 


0.343 ±0.128 


9.8 ±3.5 


66 


0.148 ±0.057 


NGC 4833* 


0.252 ±0.101 


13.6 ±5.3 


175 


0.078 ±0.031 


NGC 5024 


0.089 ± 0.040 


14.0 ±6.1 


373 


0.037 ±0.016 


NGC 5272 


0.249 ± 0.064 


41.9 ±9.9 


338 


0.124 ±0.030 


NGC 5286 


0.123 ±0.036 


23.4 ±6.6 


492 


0.048 ±0.014 


NGC 5634 


0.272 ±0.137 


14.5 ±6.4 


145 


0.100 ±0.045 


NGC 5824 


0.212 ±0.065 


41.3 ± 11.6 


529 


0.078 ± 0.022 


NGC 5904 


0.290 ±0.066 


40.1 ±8.3 


313 


0.128 ±0.027 


NGC 5927 


0.433 ±0.101 


72.0 ± 12.9 


314 


0.229 ± 0.043 


NGC 5986 


0.235 ±0.087 


36.7 ± 12.7 


403 


0.091 ±0.032 


NGC 6093 


0.152 ±0.051 


26.2 ±8.3 


381 


0.069 ± 0.022 


NGC 6101* 


0.186 ±0.097 


6.2 ±3.2 


97 


0.064 ± 0.033 


NGC 6139 


0.273 ± 0.079 


34.8 ±9.0 


302 


0.115 ±0.031 


NGC 6144* 


0.228 ±0.130 


5.7 ±3.2 


62 


0.091 ± 0.053 


NGC 6171 


0.599 ± 0.248 


22.2 ± 7.1 


65 


0.341 ±0.117 


NGC 6205 


0.218 ±0.059 


32.3 ±8.0 


436 


0.074 ±0.019 


NGC 6218 


0.312 ±0.108 


11.9 ±3.8 


91 


0.130 ±0.044 


NGC 6229 


0.273 ±0.101 


31.7 ± 10.0 


287 


0.111 ±0.035 


NGC 6235* 


0.349 ±0.163 


6.7 ±2.9 


34 


0.196 ±0.092 
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NGC 6254 


0.316 ±0.103 


17.9 ±5.0 


175 


0.102 ±0.030 


NGC 6284 


0.305 ±0.123 


21.7 ±8.2 


133 


0.163 ±0.063 


NGC 6304 


0.320 ±0.084 


44.7 ±9.8 


140 


0.319 ±0.075 


NGC 6316* 


0.272 ± 0.077 


39.1 ± 10.5 


205 


0.191 ±0.053 


NGC 6341 


0.157 ±0.062 


14.9 ±5.5 


262 


0.057 ±0.021 


NGC 6352 


0.403 ±0.172 


19.3 ±6.4 


55 


0.352 ±0.125 


NGC 6356 


0.327 ±0.063 


78.7 ± 13.6 


365 


0.216 ±0.039 


NGC 6362 


0.426 ±0.151 


18.6 ±5.5 


85 


0.218 ±0.069 


NGC 6366* 


0.398 ±0.188 


7.2 ±3.2 


20 


0.358 ±0.179 


NGC 6388 


0.290 ±0.031 


262.7 ±25.9 


1686 


0.156 ±0.016 


NGC 6397* 


0.580 ± 0.284 


7.6 ±3.5 


43 


0.177 ±0.086 


NGC 6402 


0.272 ± 0.081 


40.0 ± 10.9 


352 


0.114 ±0.032 


NGC 6426* 


0.227 ±0.117 


6.0 ±3.1 


51 


0.118 ±0.062 


NGC 6440* 


0.427 ±0.067 


123.3 ± 16.8 


423 


0.292 ± 0.042 


NGC 6441 


0.225 ±0.023 


286.8 ±27.7 


1906 


0.151 ±0.015 


NGC 6496* 


0.197 ±0.086 


9.3 ±3.9 


46 


0.201 ±0.090 


NGC 6539* 


0.495 ±0.120 


37.1 ±8.1 


117 


0.317 ±0.075 


NGC 6541 


0.237 ±0.084 


23.4 ± 7.5 


279 


0.084 ± 0.027 


NGC 6569 


0.429 ±0.109 


46.2 ± 9.9 


187 


0.247 ± 0.056 


NGC 6584 


0.214 ±0.084 


12.5 ±4.7 


125 


0.100 ±0.038 


NGC 6624 


0.433 ±0.091 


56.2 ± 10.3 


170 


0.330 ± 0.065 


NGC 6637 


0.303 ±0.080 


40.2 ±9.3 


213 


0.189 ±0.045 


NGC 6638 


0.307 ±0.101 


19.6 ±5.7 


99 


0.198 ±0.061 


NGC 6642* 


0.520 ±0.200 


11.4 ±4.0 


57 


0.201 ± 0.075 


NGC 6652 


0.416 ±0.185 


20.2 ± 7.4 


70 


0.289 ±0.111 


NGC 6656* 


0.137 ±0.054 


11.6 ±4.5 


259 


0.045 ±0.018 


NGC 6681 


0.290 ±0.103 


16.5 ±5.3 


155 


0.106 ±0.035 


NGC 6717* 


0.533 ±0.214 


9.3 ±3.5 


36 


0.258 ±0.106 


NGC 6723 


0.188 ±0.060 


18.3 ±5.4 


198 


0.092 ± 0.028 


NGC 6752 


0.430 ±0.143 


21.6 ±6.2 


181 


0.119 ±0.035 


NGC 6760 


0.470 ±0.109 


46.1 ±9.2 


159 


0.290 ± 0.063 


NGC 6809* 


0.340 ±0.131 


9.7 ±3.6 


79 


0.122 ±0.048 


NGC 6838* 


0.510 ±0.173 


13.5 ±4.3 


30 


0.452 ±0.165 


NGC 6864 


0.387 ±0.075 


66.0 ± 11.0 


369 


0.179 ±0.031 


NGC 6934 


0.211 ±0.086 


14.1 ±5.2 


194 


0.073 ± 0.027 


NGC 6981 


0.415 ±0.137 


21.3 ±6.0 


106 


0.201 ±0.060 
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NGC 7006 


0.318 ±0.071 


26.6 ±5.4 


220 


0.121 ±0.026 


NGC 7078 


0.173 ±0.056 


27.2 ±8.4 


596 


0.046 ± 0.014 


NGC 7089 


0.129 ±0.033 


35.5 ±8.7 


720 


0.049 ±0.012 


NGC 7099* 


0.147 ±0.086 


5.4 ±3.1 


130 


0.042 ± 0.024 


PYXIS00* 


0.234 ±0.172 


3.1 ±2.2 


29 


0.106 ±0.079 


RUPR106* 


0.396 ±0.218 


7.1 ±3.8 


46 


0.154 ±0.086 


TERZAN8* 


0.451 ±0.199 


7.7 ±3.2 


33 


0.233 ±0.105 
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Table 7: Other Parameters for the RGBB. N is the total number of stars used in the 
fit, B is the exponential slope of the RG luminosity function, and urgbb is the brightness 
dispersion of the RGBB in /-band. Only measurements for gold sample GCs are shown. 



Name 


AT 
iV 


7~> 

D 


<?RGBB 


LYNGA07 


329 


0.653±0.059 


0.068±0.015 


NGC 104 


2416 


0.684±0.023 


0.040±0.005 


NGC 1261 


808 


0.710±0.044 


0.027±0.009 


NGC 1851 


1241 


0.687±0.036 


0.042±0.008 


NGC 2808 


3308 


0.711±0.025 


0.092±0.013 


NGC 3201 


214 


0.594±0.073 


0.080±0.023 


NGC 362 


1060 


0.768±0.040 


0.033±0.006 


NGC 5024 


1155 


0.665±0.032 


0.028±0.011 


NGC 5272 


1325 


0.708±0.033 


0.041±0.012 


NGC 5286 


1951 


0.740±0.025 


0.028±0.007 


NGC 5634 


434 


0.716±0.057 


0.044±0.038 


NGC 5824 


1383 


0.685±0.034 


0.089±0.021 


NGC 5904 


968 


0.677±0.039 


0.037±0.009 


NGC 5927 


1103 


0.848±0.057 


0.070±0.014 


NGC 5986 


1229 


0.698±0.036 


0.048±0.025 


NGC 6093 


1286 


0.666±0.033 


0.052±0.014 


NGC 6139 


904 


0.724±0.038 


0.062±0.019 


NGC 6171 


323 


0.751±0.075 


0.101±0.035 


NGC 6205 


1252 


0.739±0.033 


0.046±0.012 


NGC 6218 


380 


0.736±0.050 


0.028±0.007 


NGC 6229 


734 


0.732±0.050 


0.061±0.028 


NGC 6254 


574 


0.783±0.054 


0.038±0.008 


NGC 6284 


644 


0.774±0.042 


0.112±0.038 


NGC 6304 


824 


0.629±0.054 


0.054±0.012 


NGC 6341 


730 


0.687±0.042 


0.030±0.011 


NGC 6352 


280 


0.688±0.105 


0.060±0.020 


NGC 6356 


1253 


0.710±0.036 


0.074±0.013 


NGC 6362 


287 


0.664±0.070 


0.060±0.017 


NGC 6388 


4003 


0.767±0.024 


0.086±0.011 


NGC 6402 


942 


0.718±0.042 


0.077±0.023 


NGC 6441 


5777 


0.758±0.021 


0.077±0.008 
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NGC 6541 


832 


0.747±0.042 


0.039±0.016 


NGC 6569 


664 


0.799±0.056 


0.086±0.017 


NGC 6584 


486 


0.741±0.055 


0.026±0.011 


NGC 6624 


892 


0.736±0.042 


0.065±0.013 


NGC 6637 


793 


0.739±0.052 


0.039±0.012 


NGC 6638 


450 


0.831±0.063 


0.044±0.012 


NGC 6652 


316 


0.842±0.101 


0.055±0.027 


NGC 6681 


448 


0.785±0.055 


0.049±0.016 


NGC 6723 


686 


0.729±0.046 


0.028±0.006 


NGC 6752 


526 


0.787±0.053 


0.056±0.017 


NGC 6760 


558 


0.689±0.053 


0.085±0.014 


NGC 6864 


1051 


0.726±0.039 


0.058±0.010 


NGC 6934 


539 


0.801±0.056 


0.029±0.017 


NGC 6981 


405 


0.674±0.060 


0.051±0.015 


NGC 7006 


761 


0.743±0.041 


0.017±0.003 


NGC 7078 


1403 


0.738±0.032 


0.059±0.017 


NGC 7089 


1855 


0.667±0.027 


0.027±0.007 



